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£52Sf2FJ2 E SIMULATION STOCHASTIQUE CENTRALISEE ET TELETRANS MISSION DE SCENARIOS 
PROBABLES POUR L'OPTIMISATION PROBABILISTE DES PARAMETRES DE SYSTEMES INDUSTRIELS 
DISTANTS. 



\gp Precede « statistique/ probabiliste » de simulation sto- 
cnastique visant & effectuer de maniere hierarchisee mais 
conjointe la prediction probabiliste du comportement futur 
et/ ou roptiniisation des parametres de rendement d'une 
multitude de system es industriels locaux distants (S u < u = 
n, >. Les systemes industriels locaux (Su) sont soumis a 

)■ 

_,, on 

genere repetitivementde maniere centralisee a des instants 
(t 0 ) successiis, une « smp serie de simulations probabilistes 
„ (W , j k n- 1..-P. j - ^r r mlj multi-valuee de description 
des variables d'etat (X , (t) <' s 1 -f» ) de I'environnement (E). 
On teletransmet repetitive ment ces « smp series de simula- 
tions probabilistes »> (Wb") en direction de Pensemble de la 
multitude des systemes locaux distants (Su). 

D 'autre part, au niveau de chacun des systemes locaux 
(Su), on collecte les « smp series de simulation probabiliste 
» (W j multi-vatuees et on effectue un traitement numeri- 
que statistique local sp£cifique sur lesdites series de simu- 
lation probabiliste (W *) recues pour predire le 
comportement futur probabiliste des systemes industriels 



locaux (Su). 

Le procede selon r invention peut etre avantageuse- 
ment utilise dans de nombreuses applications aussi bien in- 
dustrielles, financieres, qu'agricoles, etc., ou le 
comportement de sites industriels decentralises soumis a 
un meme environnement aleatoire doit etre optimise. 
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Domaone ttechmoqwe 

L'invention se rapporte au domaine technique du « Risk Management Industrie! » et plus 
sp£cifiquement : 

5 - aux processus slatistiques/probabilistes dc prediction de la loi dc probabilitc dc la 
distribution statistique des variations d'un indicateur de comportcment futur ct de 
] 'optimisation des parametres de rendemcnt d'unc multitude de systemes industriels 
locaux dtsfanis soumis a un environnement global commun aleatoire dvolutif, 

- du type par simulation, c'est4k-diie fonclionnant par dchantillonnage prealable de 
10 scenarios probables des dtats fulurs des variables d'envjronncment, 

^invention se rapporte spdeifiquement au cas : 

- oQ unc multitude de sites industriels ddcentralisds oni un comportemcnt incertain du fait 
qu'ils sont soumis a un environnement global commun aleatoire dvolutif (que Ton pcut 
dans la majoritc des cas qualifier d'environnement de risquc) > 

15 - ct ou cet environnement aleatoire dvolutif est n&mmoins quantifiable, e'est-a-dire qu'il 
pcut atre majoritairemcnt decrit num£riquement par des variables d'etat communes (au 
sens defini ci-aprds) a r ensemble des sites industriels. 
Ces variables d'etat sont en effet: 

- soil, communes a 1 'ensemble des sites industriels, 

20 - soil, lorsqu'elles sont en apparence spdeifiques a un site industricl, « quasi-commune » du 
fait que icur comportemcnt est dicte par des relations complexes avec d'aulrcs variables 
d'etat du modelc (non spdeifiques au dit site), de telle sorle que leur loi dc probabilitc nc 
pcut etre decritc independamment de celle des autrcs variables d'etat communes ou 
spccifiques d'autres sites. 

25 Par la suite, et notamment dans la description qui va suivrc ainsi que dans le texte des 
revindications, on ddsignera (par extension) du termc de « variables d'etat communes » un 
ensemble dc variables d'etat de Tun des deux types decrits ci-dessus. Par la suite, on utilisera 
dgalcmcnt mdiffercmment le terme d'environnement aleatoire ou d'environnement de risque. 
Par site indusiriel on entend, scion Pinvention : tout lieu d activity humaine ayant pour objet 

30 la production d'unc valeur ajoutde (flux de production) — ou cxprimable dircctement ou 
indirectcment en terme de valeur ajoutee — , que cclle-ci soit liee a I'industrie manufacture, 
£nergetiquc, agroalimenlaire, financifcre, de transport, de same publiquc ou d'exploitation des 
maticrcs premieres, etc. 

Par la suite on utilisera fnSqucmment commc exemplc duplication dc Invention, et comme 
35 domaine d'analyse de Tart antdrieur, le domaine financier. Car e'est dans ce domaine que 
P application des processus de prediction probabiliste de comportcmcnts aleatoires est le plus 
frequent du fail de r adoption de regie men tat ions Py contraignant (directives « GAD » du 
Conscil des Communautds Europ<Senncs, recommandations du «Comit£ dc Bale » 
rccommandalion du «Groupe des 30»): Mais les cnscignemenls de Pinvention pcuvent 
40s'apphquer dans presque tous les domaines industriels. Si Pon prend commc exemple ie 
domaine financier, on entend par risque quantifiable (environnement aleatoire ) les risques de 
marchc, de crddu, dc contrepartie, de liquidite... Par opposition, on pcut considcrcr comme 
non quamifiablcs les risques legaux, de fraudes, mediatiques, 

L'invention concernc sp^cifiqucment les cas : 
45 - ou le comportemcnt des variables d'etat de Penvironnement est incertain ou incontrolable, 
cn sorte que celles-ci constituent des « factcurs de risque » aldatoires qui pcuvent causer 
une variation et notamment une baisse d'efficacW (ou un risque sur Ic comportemcnt de 
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chacun) des sites industriels decentralises, pouvani se manifester nolamment par des penes 
(d'argent, de vie humaine, de production,...) ; 

- el oti refficacitede chaquc systfeme local peut en outre £tre decrite par un indicateur 
numerique de rendement selon un modclc numdriquc systemique local (gdndralement du 

5 type param£ire) prealablemcnt dtabli, dependant pour chacun des systfcmes, ct & chaque 
instant, de tout ou partie des variables d'etat et dc parameires Iocaux d'action permcttant 
dc regler le niveau de fonctionnement du systfcmc industricl local. 

L'invention concerne les domaincs industriels, aussi bicn manufacturiers que financiers, 

publics ou privds, od il convient : 
10 - d'estimer statistiquement pour chaque site industriel, afin de le circonscrire, un indicateur 
numerique de rendement ateatoire qui peut Strc dc type soit industriel manufacturier 
(production, nuisance), soit financier (profii/pcrtc) ou autre, et qui est Y\6 & l'excrcice 
d'une production dc valeur ajoutde (production manufacture, agricole ou financidre...), 

- et d'optimiscr Hndicatcur numdrique de rendement (ou dc risque) de chaque site 
1 5 industriel, en reglant les paramfctres locaux d'action dc cc site industriel, 

- ceci, en tenant compte notamment de entires reglcmcntaires imposanl des limitcs t eel 
indicateur de rendement (ou de risque), 

Description de Tart anterieur 

La dcuxifeme moitie du XX* mc sifcclc a vu la prise de conscience des effets pcrvers du risque et 
20 par voic de consequence de son cout el la naissancc d'une nouvcllc science: le «Risk 
Management ». Lc « Risk Management » est la science : 

- qui etudie les factcurs quantifiables ou non qui mcnacent lc rendement ou lc 
comportcment d'un systfcme industriel ou privd, 

- et qui apporle des solutions stratigiqucs (organisationnellcs ct de proeddure) ct tactiqucs 
25 (methodes quantitatives de gestion '. . .) afin de combattre ces fucleurs. 

L'art anterieur connaii des methodes de «Risk Management » de systcmes industriels. Ainsi 
le brevet US 6021402 d£crit une methode de gestion des risques d'un rescau dc systcmes 
industriels de production dlectrique. Le brevet US 5930762 decrit un systeme dc controlc de 
risque d'un dispositif dont lc comportcment est lie h des parameires multiples 
30 d'environncmcnt. 

L'art anterieur connait principalement cinq methodes diffcrcntes dc mesure ct/ou prise de 
decision relative au risque des systcmes aldatoires : 

- les processus 5 arbrc decisionnel, 

- les processus statistiques/probabilistes, (ex: Variance-covariance, Historiquc, Montc- 
35 Carlo), 

- les methodes variationnellcs (ex : La « valeur de I point de base » (PVOl)), 

- les mdthodes par scenarios extremes (ex : « Stress testing », « Scenario catastrophe »), 

- les methodes d'intelligence anificielle. 

a) Les processus h arbre ddcisionnel sont utilises dans les procedes de prise dc decision 
40 s^quentielle. Un diagramme qui ressemble aux branches d'un arbre montre Tensemblc des 
£v£nements possibles et Tenchainement des options de decision qui sont lc resultat de 
decisions pr&rddcntes. Cette mdthode h arbre ddcisionncl peut dire utilisSe uniquement dans 
les cas ires simples oto les £v£nements anticipes et la probability de chaque <5venemenl sont 
deja connus. Us calculs numeriques mis en ccuvre par cette mdthode sont simples et peuvent 
45 d'ailleurs souvent etrc implement^ sans calculateur. Mais ellc est insuffisantc pour decrire la 
structure de risque de systcmes aleatoircs complexes. 
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b) Dans dc nombreux cas, lcs evenements anticipe* dependent de difterentes variables d'etat 
donl la valeur nc peut pas £tre connue exactemcnt. Cc type de problfemes est habitueliement 
analyst par un processus statistique/probabiliste. On appellera par la suite « processus 
statistique/probabiliste » d'estimation du comportcment futur d'un systcme soumis a un 

5 environnement aieatoirc pouvant 6trc decrit a Paidc de variables d'etat, une methode 
consistant : 

- a effcctuer unc estimation statistique de la loi de probability de la distribution statislique 
des variations d'un indicateur numdrique dc rcndement dc cc systcme, 

- h partird'un module de distribution de probabilite jointc des variables d^lat de 
10 Tenvironnemcnl (factcurs dc risque) pr€alablemenl etablie. 

Ceci est en opposition avec le processus du type a arbre decisionnel (ou toute autre mcsure du 
risque) dont lcs rdsultats sbnt des valeurs uniques. 

c) Lcs methodes variationnelles quantifient la consequence sur 1c comportcment d'un syst&me 
_ dc la variation d'un ou plusieurs paramfctres d'environnement (variables d'etat). Ainsi dans le 

i 5 domaine financier, la methode dite « valeur de 1 point de base » (PVOl) quantifie la 
sensibility d'un portefeuillc dc titres financiers a un mouvement parallelc ct unilateral dc la 
courbc des taux d'intSret vers le haut ou vers le bas pour un d6ca!age de 0,01 % (un « point de 
base »). On notera que dans ce cas la probabilite d'occurrcncc et Samplitude du mouvement 
ne sont pas apprehendes. Ceci est insuffisant pour decrire la realitc du comportement de 

20 risque d*un systcme aieatoirc complexe. 

d) Lcs methodes par scenarios extremes visent a cstimcr le resultat sur un systcme dc 
J" occurrence de scenarios particuliers. Ainsi, dans le domaine financier, la methode « stress 
testing » permet dc faire varicr les taux dMntcreU a plusieurs endroits de la courbe, par 
exempic : concomitammcnt a la haussc pour une des cchdanccs et a la baisse pour une autre. 

25 On utilise cette methode a des fins de simulation de cc qui sc passerail suite a des scenarios 
catastrophes (prevision de ce qui est sense n'arriver que ires rarement). On notera que, dans ce 
cas, les aspects temps et probability n'interviennent pas. Ceci est egalemcnl insuffisant pour 
decrire la reality du comportcment dc risque d'un syslfemc aieatoirc complexe. 

e) Dans les cas differents ou Ton souhaite non pas estimcr un risque, mais plutot predire une 
30i'volution aieatoirc telle que par exemple le prix du petiole ou la variation de prix dune 

i: :ion, trop de factcurs entrcnt en jeu qui sont cux-mcmcs trop incertains el difficilcs a 
modeller avec des outils mathdmatiques. Pour aborder ces problcmcs dc prediction, on utilise 
des methodes dites «d , intelligence artificielle » a «r6seau de neurones » et des «systfcmes 
experts ». Ces techniques d'intelligence artificielle onl la capacity de raisonner sur la base 

35d'informaiions floues ou incompletes. Recemment, ces technologies ont <5t<5 utilisces dans le 
domaine de la prediction gdologique de production ct celui des predictions boursiercs. Mais 
elles sont inadaptees pour des applications generates d'analysc de risque de systcmcs 
tndustricls, cu 6gard a des critcres (notamment rdglcmcntaires) communs. Ainsi le brevet US 
5521813 decrit un systfcme de prediction des consequences des variations climatiques sur un 

40 systcme produclif, utilisant une methode a regression multiple et rcscau de neurones. 

La presente invention se rapporte specifiquemcnt au domaine technologiquc des processus dc 
statistiqucs/probabilistes qui visent a effcctuer une estimation statistique de la loi dc 
probabilite de la distribution des variations d'un indicateur numeriquc de rendement de ce 
systcme, a partird^un modele de distribution de probabilite jointe des variables d'etat de 

45 renvironncmcnt. Le domaine de Tart anterieur le plus courant dans lcquel lcs processus 
statistiques/probabilistes d'analyse de risque sont utilises est eclui de la finance ct plus 
specialement eclui des risques de marche. Dans eel environnement global de Pindustrte 
flnancfcre, le risque de marchd est le risque auquel est confronte Tinvcstisseur en raison de sa 
meconnaissancc de revolution future des variables d'etat communes fondamcntales du 

50 marchd que sont : les cours des titres, les taux d'interets, les taux dc change.... Ces variables 
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de marchd, ygalement appe!6es facteurs de risque, determinem le prix des litres, des aciifs 
contingents, dcs portefeuilles, ...etc. 

Lcs risqucs associ^s & differents produits financiers pcuvent etre lvalues par quelques 
nombres caractiristiques : duration et convexite pour lcs obligations, beta et performance 
5 espdrde pour les actions, valeur-temps, volatility, delta, gamma, theta, v^ga (ddrivgc par 
rapport a la volatilitc), etc., pour les options. Mais ccs diffcrent.es approches du risque dites 
par « sensibilite » : 

- sont specifiques h un type de produit financier, 

- s'appuicnt sur des hypotheses de distribution specifiques, 
10 - ne sont pas faciles h interpreter, 

- sous-estiment souvent le risque global de pcrtc sur un portcfcuillc qui est pourtant celui 
qui importe pour Tinvcstisseur, 

- n'^valuent pas le niveau de probabilite de survenance de ccs risqucs, 

- et ne prennent pas T horizon temporel en consideration. 

15 II est preferable d'avoir unc approchc unificc probabiliste et temporcllc du risque. Ainsi, Tart 
anterieur utilise la methode dile « Value at Risk » (VaR) commc indicateur de risque d'un 
systeme financier soumis & un environnement aldatoirc de marchy. Pour un actif financier (ou 
un portefeuille financier), la VaR est un montanl tel que la perte encourue sur cet actif (ou ce 
portefeuille) durant un certain intervalJe de temps ne depassera VaR qu'avcc une probability 

20 de niveau (faible) donne. Le calcul de la VaR necessite la connaissance de la loi dc probability 
de la perte subie (indicateur numerique dc rendement). Le probleme fondamental des 
processus stalistiques/probabilistes d 1 analyse dc risque (et notarnmcnt dcs processus de calcul 
dc la VaR) est done de prevoir la distribution de loi de probability de la distribution statistiquc 
des variations d^un indicateur nurn<5rique de rendement (pcrtc dans le cas de la VaR). 

25 Lcs processus statistiques/probabilisles d'estimation statistiquc de loi de probabilite dc la 
distribution dcs variations d'un indicateur numerique de rendement d'un sysifcme son! 
principalement de deux types : 

- Ccux, que nous appcllcrons par la suite « par simulation », 

• qui procedent par gdndration d'une multitude d'echantillons statistiques multi-varies 
30 ( dc i-uplcts) de description de valeurs probables conjointes dcs variables d'ytat 

(parametres de risques), 

• qui abouttssent h la description « de Thistogrammc point par point » de la structure dc 
la loi dc probability de la distribution staiistique dcs variations de Tindicatcur 
numerique dc rendement du systeme, (e'est-a-dire dans son integrality), 

.15 ♦ et qui pcrmettent notamment de calculcr la VaR de manure gcncralc, sans hypothese 
prdalublc sur la loi dc probability de eelte distribution statistiquc dcs variations de 
Tindicatcur numerique de rendement ; 

- Par opposition 5 ceux, que nous appellerons par la suite « analytiques » (et que Ton 
pourrait egalemem qualifier dc parametriques, normatifs ou synthytiqucs,. . .)> 

40 • qui ne fournissent que quelques parametres (cn nombre fixy) descriptifs de la loi de 
probabilite. de la distribution statistiquc des variations de T indicateur mimyrique de 
rendement du systeme, 

• qui imposcnt unc forme particuliere au module devaluation de Tindicatcur numerique 
de rendement du systeme en fonction dcs variables d'ytat (paramfctres de risques), 

15 (notamment une hypothese de lindarity, quadralicile, etc.), 

• et qui aboutissent a des hypotheses iimitatives dc description parametrique spycifique 
de la loi de probabilite de la distribution statistiquc des variations de Tindicateur 
numerique de rendement. 
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Ainsi, Tart anteneur connaii principalement qualrc mtflhodes classiqucs d'cstimation de la 
VaR (et plus gdneralernent d'estimation statistique de loi de probability de la distribution 
statistique des variations d'un indicateur numdriquc dc rendement dun systcme soumis h un 
environnement aldatoire) : 
5 - deux m&hodes statistiques/probabilistes analytiqucs : 

• la technique « delta normal » ou technique « de la matrice des van ance-co variance 
estim£es », qui suppose le rendement 1in£aire ; 

• la technique « delta-gamma », qui tient compte d'un terme quadratiquc dans le calcul 
du rendement du systcme ; 

1 0 - et deux mcthodes de type statistique/probabiliste par simulation : 

• la technique de simulation de Monte-Carlo ; 

• la technique de r analyse historique. 

L' invention sc rappoite exclusivement aux domaincs techniques des procedds de type 
« siatisliques / probabilistes par simulation ». 

15 Pour le calcul de la VaR d'un portefeuille financier, les entries sont classiquement : 

- les modeles mathematiqucs devaluation dc prix des diffcrcnls actifs, ex : 

• 1c modele de prix par combinaison lineairc de variation du prix d'un actif en fonction 
des variations des facteurs de risque (tel que le « module h indice simple de Sharpc » 
pour les actions, qui postule que la variation du prix d'une action est une fonction du 

20 premier degrd dc la variation d'un indice gcndral du marche*), 

• les modeles de prix non lin<5aire, tel que le module en « dcvcloppemcnt dc Taylor » 
independant du temps, 

• et les modeles dc prix fonction du temps, notamment du type modele aleatoire 
stationnairc, tel que « ]e modele de Black ct Sholcs » pour revaluation des options ; 

25 - et les hisloriques, c'esi-a-dire les observations des prix des differcnts actifs pour un 
certain nombre de periodes passtes. 

La methodc « delta normal » (ou linemre) ditc de la matrice des varianccs-covariances 
estimecs est cell e qui a ele proposed ct ddvcloppcc par JP Morgan, avee son systeme 
RiskMetnx . EIlc consiste en trois 6tapcs : 
30- le relevd des facteurs de risque elcmentaircs en lesquels les actifs financiers du 
portefeuille pcuvent etre decomposes, 

- la repartition des flux financiers assoctes h ces facteurs de risque elcmentaircs en fiux 
agregds correspondant a des defiances standards, 

- 1c calcul effectif de la VaR. 

35 Cette methodc. fait les hypotheses theoriques specifiqucs limitatives suivantes : 

- les facteurs de risque lids h chaque actif obcisscm & une loi normale, 

- la performance des actifs depend lineairement des facteurs dc risque, 

- ct par voic dc consdquence la performance globale du portefeuille depend Iiniairement 
des facteurs de risque. 

40 Le calcul de la VaR s'effectuc alors a partir de la matrice M dite de « variance-covariance », 
de dimension nxn dans le cas de n facteurs de risque, et constitute & partir des /ix(/i+l)/2 
variances el covariances des variations des facteurs dc risque, calculees sur la base d'un 
historique passe\ Si Ton appelle (X = A/'* '■ •■">) les scnsibilites du portefeuille sur les different* 
facteurs de risque (X = Xr l *" > ), la VaR est donncc par la formulc : 

oh z est le quantile de la distribution normale rcduitc pour la valeur de rintervaHe de 
conuance de la VaR rctenue. En sorte que, selon la methodc dc la matrice de variance- 
covanancc, on calcule la structure du risque, dans un environnement h n facteurs de risques, 
en gdndrant un ensemble de iT paramctres fixes dc la matrice de variance-covariance 
50descnptifs des variances ct covariances des variations des facteurs dc risque. C'est en ce sens 
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que la methode de la malrice des variances-covarianccs est une m&hode 
statistique/probabilistc analytique. L'hypothese de lindaritc retenue necessite pour etre mise 
en oeuvre des hypotheses strides sur les sensibilites du portefcuille, la rendanl inapplicable 
pour des portcfeuilles incluant des risques dits « gamma » ou dc « convexite ». 

5 La methode dite « delta-gamma » est une forme plus dlaborce de la methode de la matrice des 
variance-covariarice. Elle abandonne l'hypothese de lincaritc an profit d'unc hypothcse 
quadratique. Ellc incorpore un second ordre dc sensibilite dans les analyses. Elle n6cessitc la 
mise en ceuvre de calculs malhematiques sophistiqucs. Malgre sa superiorite vis-a-vis de la 
methode « delta normal*, la methode « delta-gamma » est inadaptde dans le cas general 
lOnotamment pour evaluer le risque de portefeuillcs incluant des options arrivant a maturite ou 
des derives exotiqucs tcls que des options « knock-out ». 

Au contraire, scion la methode de simulation de Monlc-Carlo. prince notamment par Bankers 
Trust, avec son systfeme RaRoc 2020™: 

- On fait le choix d'une famille de modeles devaluation pour les differents facteurs de 
1 5 risque ou les actifs du portefeuille, justifies a partir d'obscrvations historiqucs. 

- A partir des historiques des actifs en presence et des facteurs de risque, on estime : 

o la distribution des differents facteurs dc risque, ainsi que des parametres qui lui sont 

assocics : moyennes et malrice de variancc-covariancc . . ., 
o ainsi que les relations liant les prix des actifs et les facteurs de risque. 
20- Pour Penscmble des differents facteurs de risque et sur la base des modeles probabilistes 
retenus el des parametres de distribution determine*, on emet grSce h un gendratcur 
pseudo-alcatoire, un grand nombre (quelques centaines a quelqucs millicrs) d'echantillons 
pseudo-aleatoircs extraits d'une distribution estimcc des variations probables futures des 
facteurs de risque. 

25 - On introduit les valeurs des echantillons de chacune de ccs simulations dans les relations 
dc comportcment des actifs en fonction des facteurs communs de risque. 

- On reconstituc pour chacun des echantillons dc simulation la valcur correspondante du 
portefcuille. 

- On construit Phistogramme global de loi dc probability de la distribution statistiquc des 
30 variations de valeur du portefeuille, sans hypothcse prealable sur cetlc derniere. 

C est en ce sens que la methode dc simulation dc Monte-Carlo est unc methode 
statistique/probabiliste par simulation. 

La m&hodcd'analysc historique, prftnte notamment par Chase Manhattan, avec les systfcmes 
C harisma et Risk$ est une simplification de la methode Monte-Carlo. Selon cetle 
35 methode d analyse historique, on emet la multitude d'echantillons des variations probables 
tuturcs des facteurs dc risque dircctement a partir dc Phistoriquc du marchc. A litre 
d exemplc, la VaR historique h un jour d'un ponefeuille pcut etre estimee a partir de cent 
echantillons de scenarios de la mani6re suivante : 

- On collecle les valeurs reconstitudes des facteurs dc risque h partir de Phistoriquc du 
40 marche sur les cent dcrniers jours de eolation. 

- Pour chaque jour, on determine le pourcentagc de variation de chacun des facteurs de 
risque. Ccla permct de simuler cent echantillons multi-values de valeurs probables de 
chacun des facteurs de risque. 

- On introduit les valeurs des echantillons corrcspondant it chacune dc ces observations 
45 historiqucs dans les relations de compoitemcnl des actifs en function des facteurs 

communs de risque. 

- On reconstitue pour chacun des echantillons historiqucs la valeur correspondante du 
portefeuille. 

- On construit Phistogramme global de loi de probability de la distribution statistiqne des 
50 variations de valeur du portefeuille, sans hypothcse prealable sur cettc derniere. 
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C'cst en ce sens que la methode par analyse historiquc est dgalement unc mfthode 
statistique/probabiliste par simulation. 

Lc principe de Tanalyse de risque utilisant unc mclhode << statistique/probabiliste par 
simulation » est connu de Tart anterieur et utilisd dans de nombreux domaines aussi bien 
5 financiers qu'industricls. La presente invention qui comporte (en combinaison) la mise en 
oeuvrc d'une mdthode de ce type ne revendique nullcmcnt ce moyen pris seul. On peul citer le 
brevet FR2769703 qui decrit un proeddd pour determiner scion la methode Monte-Carlo 
Pincertitude dc mesure d'unc machine & mcsurer en coordonndes. Le brevet EP0837408 decrit 
un proedde de simulation de pulverisation par mclhode Monte-Carlo. Le brevet W09905685 

10 decrit un systeme automatique de surveillance industriclle par mdthode Monte-Carlo. Le 
brevet W09859322 decrit un systfcme dc gdndration des valeurs de pixels d'une image par 
mdthode Monte-Carlo. Lc brevet W09852132 decrit un systeme d^affectation dc rcssources 
organisationnelles el devaluation des couts et prix des ressourccs par methode Monte-Carlo. 
Le brevet WO9850835 decrit un controleur programmable pour supcrviser les activity 

15 d'dquipemenls industriels par mdthode Monte-Carlo. Le brevet W09836364 ddcril un 
processus devaluation d'intcgrale par methode Monte-Carlo. Le brevet US5940810 decrit 
une methode d'estimation dc la valeur d'actifs financiers complexes par mdthode Monte- 
Carlo. 

On constate que Ton pcut rdpartir les mdthodes statistiqucs/probabilistcs d'estimalion de loi 
20 de probability dc la distribution statistique des variations d'un indicateur numdrique de 
rendement d'un systeme soumis h un environnement aleatoire (ct notammcnt les methodes 
d'estimalion de la VaR) en deux groupes distincts: 

- les methodes analytiqucs selon lesquelles on gdndre un nombrc fixe ct petit dc paramctres 
descriptifs du comportemcnt statistique des param&trcs de risque, ct aboutissant h unc 

25 description paramdtrique dc la loi de probability dc la distribution statistique des variations 
de Findicateur numdrique dc rendement d'un systeme, (e'est le cas notammcnt de la 
mclhode de la matrice des variances-covarianccs dc RiskMctrics ™ « delta normal » ou 
« delta-gamma ») ; 

- et les methodes par simulation qui procfedent par « dchantillonnage de scenarios probables 
30 de variables d'envjronnemcnt » et selon lesquelles on genfcre unc multitude (donl le 

nombrc peut ctre choisi) de tirages repetds d'echaniillons dc valeurs futures probables 
plausibles des variables d'environncment, aboutissant h unc description globale absolue de 
rhislogramme de loi dc probabilite de la distribution statistique des variations d'un 
indicateur numdrique de rendement d'un systeme (e'est le cas notammcnt de la methode 
35 de Monte-Carlo ct de la mdthode d'anaiyse historique). 

L'invention sc rapporle exclusivemcnt aux procedds d'analysc de comportement : 

- de type « statistique/probabiliste par simulation », 

- effectuant un « dchantillonnage (tirage) de scenarios probables multi-valuds des variables 
d' environnement ». 

40 1/ art anterieur sail, au niveau d'un site unique, metlre en ceuvre une mdthode de type 
« statislique/probabiliste par simulation », et notammcnt effectuer Tcnscmblc des acquisitions 
dc donndes, calculs statisliques, etc... en vue dc determiner la loi de probability dc la 
distribution statistique des variations d'un indicateur numdrique de rendement d'un systeme 
lorsqu'il est soumis 5 un environnement aleatoire. L'art antdrieur sait dgalcment calculer la 

45 VaR du portefeuillc d'unc institution financterc seule, par mdthode 
« statislique/probabiliste par simulation ». 

n est par ailleurs connu par Tart antdrieur d'dmettrc rcpetitivemcnt des valeurs numeriqucs 
d'un site centralisd en direction d'une multitude de sites distants. Cest notamment le cas des 
ecrans de donndes financidres de type Reuters ou Bloomberg. Un site central est relic h des 
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sources de donndcs variables ei aleatoires d'environnements (cours d'actifs financiers, etc.). 
II concentre ces donndcs et les 6met rep£tilivcment dc mantere decentralisce cn direction 
d'une multitude dc sites distants (ecrans d'opdratcurs financiers), Mais on constatera que 
selon Part antericur, les donnees t^letransmises sont dcs donn£es rdelles observables (valeurs 
5 de titres) ou sont quasi-observables, c'est-a-dirc reconstitutes par calcul (valeurs d'indiccs, 
taux de rendement internes...) et destinies h etre utilisdes comme telles. L'art anterieur ne 
connait pas 1c principe de la ^transmission systematique et repetitive d'dchanlillons 
probables multi-valuls de valeurs futures possibles de variables d'environnemcnt cn vuc de 
reconstituer un hisiogramme. 

10 Les methodes de Tart anterieur telles qu'exposecs ci-dessus nc pcrmcUcnl pas d'effectuer de 
manifcre satisfaisante la prediction probabilisie du comportement futur elA>u r optimisation 
des paramfctres de rendement d'une multitude dc systfcmes industriels locaux distants soumis 
a un environnement global commun. Notammcnt elles ne pcrmcttent pas de calculer 
concomitamment dans de bonnes conditions la VaR d'une multitude dc sites financiers 

1 5 dislincts gerant des ponefeuilles dc titres. En effet : 

D'une part, comme on Fa vu plus haut, seuls les processus statistiqucs/probabilistes qui visent 
a effectucr une estimation statistique de loi dc probability de la distribution statistiquc des 
variations d'un indicateur numerique de comportement, donnenl une mesure satisfaisante 
fiablc et fidele des risques. 

20 D'autre part* parmi les processus slalisttques/probabilistes : 

- A) La methode analytique, et noiamment la methode dc la matrice de variance-covariance 
estimae, est simple el necessite peu de calcul. Ellc pcul d'ailleurs etre mise en oeuvre 
conjointemcnt au niveau dc plusieurs sites distants. Ainsi, selon une variance dc mise en 
ocuvrc dc Part ant£rieur, 1c systeme RiskMetrics ™ fournit a partir d'un site centralise, par 

25 Internet, et en direction d'une multitude de sites distants qui utiliscnt sa methode, unc 
emission repctitives d'un nombre fixe de parametres de la matrice dc variance-covariance. 
En sorte que le calcul global dc la VaR est hicrarchisc cntrc 1c site central et les sites 
distants. Par contre les valeurs transmises ne sont pas dcs echanliilons probables multi- 
values de variables d'environnement mais de simples parametres (voir plus haut). Le 

30 defaul de cette methode est quelle necessite des hypotheses dc base trfcs strides 
(normality conditionncllc, linearite, stationnarite) qui nc sont pas satisfaites pour une 
grandc partie des instruments financiers. En sortc que lc calcul de la VaR selon la mdthode 
de la matrice de variance-covariance produit dcs erreurs dans Ic cas de ponefeuilles 
comprcnant des instalments ayant un « gamma*, unc « convexity » (voire des inversions 

35 dc convexity ou des discontinues) tels que les options, les notes structurees, les 
« mortgage-backed securities*.,. Plus gdneralement, les VaR analytiqucs, y compris les 
plus pcrfcctionndcs, rdsument le comportement d'une structure dc risque a remission d'un 
nombre rcstreint de paramfetres, ses « grecques » (delta, gamma, vdga, theta, ...). Ceci est 
cn general insuffisant pour la description d'environnements ayant dcs facteurs dc risque 

40 complexes. 

- B) La methode classiquc par simulation, et notammcnt la technique de simulation dc 
Monte-Carlo est consideree par la majority des specialises comme la plus « pure », la plus 
globalc, « la plus toum£e vers Tavenir » dans sa conception. Elle est incontestablement la 
plus juste theoriquement, la plus precise et elle s'applique h la totalitc dcs instruments 

45 financiers. Cepcndant cllc necessite dcs choix complexes de modeics devaluation ainsi 
que des processus stochastiques de gyration d'echanlillons aleatoires correctcmcnt 
choisis. Ceci est tres complcxe €tant donne la difficult de revaluation de certains 
instruments ct leurs constantes Evolutions, Elle mobilise des equipes de spccialistcs 
nombreuses pour etre correctement impldmentce sur un site. En outre elle necessite pour 

50 etre precise la gdndration d'un nombre tres important de tirages pseudo-alcatoires 
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d'eehantillons. Ceci rend necessaire la mise a disposition dc sysicmcs informatiques de 
grande taille et coGteux. En resume, scion la mcthode Monte-Carlo de Part anterieur, 
chaque site doit disposer d'un systeme informatique couteux et d'une large et cofiteuse 
equipe de mise en ceuvre. En sorte quelle est pcu accessible a la majority des institutions 
5 financiers ou sites industriels qui pourtant souhaiteraicnl cl auraicnt int6r£l a la mettre en 
ceuvre pour rdduire leur risque. 
- C) LMnconvdnicnt majcur dc la mcthode dc simulation par analyse historiquc est sa tres 
grande sensibility a la qualite des donnees. En effel, une ou quclqucs observations sonant 
dc Pordinairc (qu'elles soient vrairnent cxccptionncllcs ou qu'unc crreur soil a leur 

10 origine) influencent fortement la valeur de la VaR pendant une longue periode (egale a la 
dur6e des historiques). On a dit que la theorie des valeurs extremes pouvait pallier cet 
inconvenient. Malheureuscment, la masse des calculs que son application entratne est 
prohibitive. En outre, les VaR historiques sont imprcciscs et doivent 6tre completes par 
des « stress tests » adequats. Par consequent, selon la mdthodc historiquc dc Part anterieur, 

1 5 chaque site doit disposer d'un systeme informatique cofltcux cl d'unc large equipe dc mise 
en ceuvre. En sorte qu'elle est dgalemcnt pcu accessible a Iu majority des institutions 
financieres ou sites industriels qui pourtant souhaiteraient et auraient int<5r6l a la mettre en 
ceuvre pour reduire leur risque. 

II est connu par Part anterieur de diffuser un scenario le plus probable dc variables 
20d , environnement en direction de sites distant* (cx : lelcviscurs). C'cst le cas notamment des 
previsions de variables metcorologiques (vents, temperatures ...) en plusieurs points. Mais 
ceci nc pcrmet pas devoir une vision probabiliste globale (loi de probabilitc de la distribution 
statistique des variations d'un indicateur numdrique dc comportement) mais limitee h un petit 
nombrc dc paramctres au niveau de chacun des sites, 

25 En oulre, les procedures de contrite des risques rcquicrcnl unc uniformisation des mesures de 
risque, done une unite de methode dans leur evaluation. 

En sorte que Panalysc dc Part anterieur fait ressortir le besoin d*unc technique dc prediction 
probabiliste du comportement futur (et/ou d'optimisation des paramctres de rendement) 
d'une multitude de systemes industriels locaux distants (notamment d'unc multitude dc sites 

30 financiers) soumis a un environnement global commun, donnani une mesure satisfaisante 
fiable et fidclc des risques, de type « statist ique/probabilistc par simulation » e'est-a-dire 
fonciionnant par dchantillonnage de scenarios probables de variables d Environnement, en 
utilisant globalement le minimum : 
de moyens humains, 

35 - et de moyens en materiel de calcul informatique. 

Resume rile Poowention 

Lc Proccdc « statistiquc/probabilistc » dc simulation stochastique selon Pinvcntion vise a 
effectuer de manicrc hierarchisee mais conjoinle (e'est-fi-dire mutualisdc) la prediction 
probabiliste du comportement futur et/ou Poplimisation des parametres dc rendement d'une 
40 multitude de systemes industriels locaux distants (S u (u = '■■ n \ II s'applique au cas : 

- ou ces systemes industriels locaux (S u ) sont tres distants les uns des autres, mais sont 
soumis a un environnement global commun (E) evolutif, pouvant etre decrit 
numenqucment par p variables d'etat (*,</)" = 1 ■ consliluanl les composantes d'un 
vecteur d'etat (X(/) = Xfyf = y ' p \ 

45 - ou les p variables d^lat (Xftf* sont communes a la multitude des // (n etant tres 
grand » 1) systemes industriels locaux distants (S H lus u ' n) ) decentralises, 

- et ou Pefficacite de chaque systeme local (5 w (wa 1 n) ) peut en outre etre decrite par un 
indicateur numcrique de rendement (lj as] - "\ selon un modcle numerique systemique 
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local paramdtrd (D u ) prdalablement etabli (/ M = D M (Mssl ~ M) X(f))), dependant pour 
chacun des systemes (S„) et ^ chaque instant (/) : 

o de tout ou partie dcs variables d'dtat (XM (im 1 " p) ) du vecteur d'etat (X(/)), 
o et d'un vecteur de^y parametres (reglables) locaux d'action (\ u = (A,,,, ^)). 

5 Le procddd dc simulation stochastique selon Tinvention procide selon unc mcthode classique 
de type « statistique/probabiliste par simulation » qui consistetout d'abord & effcciuer un 
«dchantilIonnage de scenarios probables des variables de Tenvironncment », puis h effecuier 
une estimation statistique de la loi de probabilitd dc la distribution des variations d'un 
indicateur numdriquc de rendement de chacun des systemes industricis locaux (S t( ) h pariir de 
lOcet echantillonnage de description de la distribution dc probabilitd jointe des variables d*6tat 
de Penvironncment prdalablement dtabli, 

Mais 1c procedd selon I 'invention vise en outre : 

- h centraliser (mutualiser) et a mettre en commun au sein d'un site central les moyens 
lourds (cn materiel informatique et en personnel) nccessaircs & la gdndration de ce( 

1 5 dchantillonnage de scenarios probables des variables de I'environncmen!, 

- et & ne realises au niveau de la multitude des systemes industriels locaux, que des calculs 
simples et peu ondreux en materiel informatique et cn personnel specialist 

Dans sa forme la plus gdndrale, le procedd dc simulation stochastique selon V invention est 
rcmarquable par la combinaison des deux moyens suivants: 
20 □ D'une part, au niveau d'un site centralist dmctlcur (O), trcs distant (de plusicurs 

kilometres, voire dc plusieurs centaines dc kilomfctres) d'une majorite des n systemes 

industriels locaux distants (&/"" 1 ,n) ), 

0 on gdnfcre repdlitivement de manierc ccntralisee a des instants (/ 0 ) successifs, unc 
«smp scrie de simulations probabilistes » (W,/ u=l ■ ~*< ks] --*>) mulli-valuee de 
25 description des variables d'etat (*/(/)'' = , ' '' ) ) dc Tcnvironncment (£). Ccttc « smp 

sdrie dc simulations probabilistes » est, scion I' invention, constitute : 
o soil d'une sdrie de « pseudo-dtats probables futurs » dcs variables (X t itf" 1 " p) ) 

de Tenvironnement (£*), 
o soil d'unc sdrie dc « pseudo-variations probables futures » dc ces variables 
30 (AX,iO t,ml '" p \ 

(ces sdries dtant simulces k un (ou plusicurs avec m > I) horizons dc temps 
(,,«-»•■■»> >fo); 

o et on tdldtransmet rdpdtitivement ces « smp series de simulations 
probabilistes »(W i J i ) multi-valude dc description des pseudo-dtats . ou pseudo- 
35 variations probables futur(e)s des variables de Tcnvironnemcnt, ce de mantere 

rcpartic ddcentralisde : 
o du site emettcur(L0, 

o en direction de Pensemble de la multitude dcs systemes locaux distants (S u ). 
a D'autre part, de mantere independante, au niveau dc chacun des systemes locaux (SJ, 
40 majoritairemcnt tres distants du site centralise dmelleur (£/), rdpetitivement : 

o on collecte lesdites «smp series de simulation probabiliste » (W,/) multi-valudes de 

description dcs pseudo-dtats ou pseudo-variations probables futur(e)s, 
o et on effeclue un traitcment numerique statistique local specifique sur lesdites sdries 
de simulation probabiliste (IV,-/) revues pour prddire lc comportcment futur 
45 probabiliste des systemes industriels locaux (5„). 

Eventuellemcnt, on rdglc lc niveau du vecteur dcs parametres locaux d 'action (k u ) = (A„|, 
...,A,, f/ ) de chaque systeme local (S u ) & une valeur optimalc calculdc en fonction dcs criteres 
prcctablis, afin notamment d'anticiper le comportemcnt probable futur cl Tefficacite de 
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chaque systfcme (S u ) soumis h Penvironnemcnt (E) au cours du ou des horizon(s) de temps (/*) 
retenu(s). 



Figures 

Sur )es dcssins : 

5 - La figure I : d6crit une rruSthode classiquc rccommandee par Pinvcntion pour la saisie de 
mcsures observables de Peiat d'un environnement et leur transmission par un fournisseur 
de ces donndes observables h un site Smetteur dc donndes probabilistes selon rinvenlion. 

- La figure 2 : ddcrit les principals etapes du processus de gyration d'un dchanlillonnagc 
dc scenarios probables dc variables d'etat d'environncmenl par un site emetteur de 

10 donates probabilistes selon l'invention procedanl par mdthodc de Monte Carlo. 

- La figure 3 : decrit les principals dtapes du processus selon Pinvcntion dc iraitemcnt de 
sdries d'echantillonnages de scenarios probables de variables d'etat d'environncment 
rc^us par un groupe de trois syst&mes industriels locaux distants fonclionnant selon 
Pinvcntion, Tun des trois systemes etant ddcrit en detail. 

15- La figure 4 : decril globalement les principalcs dtapes du procede selon r invention de 
simulation stochastique repetitive hierarchisec par echantillonnage de scenarios probables, 
pour la prediction probabiliste du comportement futur ct/ou Poptimisalion des parametres 
de rendement d'une multitude de systdmes industriels locaux distants. 

- La figure 5 : decrit globalement les principalcs Stapes d'une variance du proedde scion 
20 Finventionde simulation stochastique rdpdtitive hiSrarchiscc par echantillonnage de 

scenarios variationnels probables, avec tcleadaptation des parametres num6riques de 
calcul du rendement dc chaque systeme industriel distant. 

- Les figures 6, 7 et 8 : decrivent trois methodes pour effectuer selon Pinvcntion la 
transmission dc «smp sdries probabilistes » du site centralise <5mctlcur vers la multitude 

25 des sites industriels locaux distants et trois m&hodes pour filtrer sclcctivcmcnt le contenu 
des « snip series probabilistes » transmises. 

- La figure 9 : ddcril une application parlieultere du procede selon rinvenlion dans le 
domaine de r Industrie financi£re. 

- La figure 10 : decrit une variante de mise en ccuvre de Pinvcntion pour effectuer une 
30 optimisation probabiliste groupie de certains groupes dc systemes industriels locaux. 

Description de reinvention 

Les figures I, 2 et 3 illustrent Porganisation gdndrale des principalcs Stapes du proced<5 par 
simulation stochastique de generation r£p6titive ccmralisce d'echantillons statistiques et de 
tdldtransmission dc scenarios probables, pour hi prediction probabiliste du comporlemenl 

35 futur el/ou Poptimisation des parametres de rendement d'unc multitude dc systdmes 
industriels locaux distants (S u ). Dans Pcxemplc rctenu, les systfemcs industriels locaux 

distants (S u ) sont consumes par des centrales electriques (5 ,&,) places dans un mcme 

environnement (£) et formant, dans ce cas particular, un rescau dc distribution clcctriquc 
d'un etat ou d'une region. 

40 On voiC.cn reference aux figures 1 et 4, Penvironnemcnt (£) d'usines de production 
electrique, connectde & un meme rdseau, influencees par des phenom&nes aleatoires que sont 
le climat (1) agissanl directcmcnt sur les performances des usincs (rendement des centrales, 
niveau des barrages), les caractcristiqucs dc la demande (quantitc, prix) (2) affectant les 
revenus de la centrale electrique, la disponibilite du nSscau (3) influemjani ie regime de 

45 fonctionnement des usines, et les caractdristiques des Energies primaires (4) influcn^ant les 
couts d'approvisionnement. Ces phenomenes aleatoires lies au comportement de 
Penvironnemem (E) du r&scau sont donnas a titre d'exemple. Mais en comprendra que dans la 
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pratique les phdnomfcncs aldatoircs h prendre en compte sont beaucoup plus nombreux et 
peuvent etre diffcrcnts. 

Meme lorsqu'ils sont locaux (commc les phdnomfcnes clima(iques), ces ph<5nomfcncs 
alcatoires nc peuvent etre considers inddpendamment des phenom&nes communs & loutes les 
5 usines (comme le prix du pdtrole, ]e niveau de la demande sur le rescau). Ainsi, la temperature 
h proxiniitc d'une usine agit non seulement sur le rcndcrncnt de F usine (done sur son cout 
marginal), mais interagit avec le niveau dc la demande et le prix de FdlcctriciuS. On a de fortes 
chances dc voir evoluer celle temperature locale de manterc corrclce avec la temperature h 
proximitc des autres usines, et done de voir le cout marginal dc ces dcrnifcres varier dans 1c 
10 meme sens et interagir par aillcurs avec le niveau de la demande. Rechercher 5 optimiscr la 
marge d'unc usine en faisant une analyse siatistiquc des parametres locaux de climat 
ind<5pendamment de r analyse -statistique sur les paramfctrcs globaux conduirait h des erreurs 
manifestes d'oplimisation. Car parametres locaux et parametres globaux interagissent dans le 
cadre d'un meme environnement global commun. 

15 En reference aux figures 1, 4 et 5, on constate qu'un centre fournisseur de donndes (Fo) est 
charge de reunir des donnees observables pennedant de decrirc F£tat de Fenvironnement 
global (E) du r&seau electriquc. Dans ce cas particulicr dc F Industrie eleclrique, le centre 
fournisseur de donnees (Fo) peut faire partie de la compagnic gdrant V organ i sat i on du reseau 
eleclrique. Mais dans d'autrcs applications, notamment financifcres, ce centre fournisseur dc 

20 donnees (Fo) est une organisation independantc (telle que les compagnies Reuters, 
Bloomberg) fournissant commcrcialement des donndes observables it une multitude de clients. 

On a installc en differents points dc Fenvironnemenl (E), une s^rie de r captcurs (y, t <h=] '" r) ). 
lis constituent des observateurs de F6tat d'environncmcnt (E). Le fournisseur dc donnees (Fo) 
collccle h distance, h chaque instant (/), les r mcsurcs de Fetal de Fenvironnement 
25 (///,(/)' * ~ r) ) en provenance des capteurs. Ces /• mesures constituent, h chaque instant 
d'tfchuntillonnagc (/*), les composantcs d'un vecteur tcmporel de mesure d'cnvironncmcnl 
(M(/*)) = (/i*(/*) tt " , "" r> ). Le fournisseur de donnees (Fo) est <5quipe d'un gmcttcur (Ao) dc 
donnees observables. Celui~ci dmct sur un reseau (Ro) demission de donnees observables : 

- soit directement les mesures recucillics par les capteurs (variables observables), 

30 - soil les rcsullats de calculs alg^briques d&erministes (tcls que moyennes, ...) op<5r<5s par le 
fournisseur dc donnees (Fo) sur ces mesures (variables quasi-obscrvabJcs). 
Par la suite, on assimilera variables observables et quasi-observablcs que Ton appcllera 
globalemenl « variables observables ». Dans le cas represent figure 1, le r<5scau (Ro) 
d'emission de donnees observables est schematise par une antenne parabolique demission (5) 

35 transmcttanl les donnees (mesures) observables (M(r*)) = (il*(/*) f *" 1 r) ) de renvironnement 
(F) h un satellite (6). 

Un site centralise emetleur (U) selon Finvcntion est decrit dans sa forme gdnerale en 
reference h la figure 2. II est generalemcnt constitu<* d % un systeme informatique (7) 
comprcnant principalcment : 

40 - une interface dc mod£Iisation (CI) permettant & un operatcur matWmaticien de faire le 
choix des diffcrcnts modeles et variables, 

- un calculateur (C) de simulation pour effectuer les calculs num£riqucs et statistiqucs, 

- une memoirede masse (#) centrale de donnees probabilistcs, 

- un rdcepteur (Bo) dc mesures observables, 

45 - et un <Srnetteur (Ap) de series dc simulations probabilistes. 

Ce site centralise emetteur (U) est tr£s distant (dc plusieurs kilometres, voire de plusicurs 
centaines de kilometres) d'unc majorite des n systemes industriels locaux distants (5,/ M= r n \ 

les centrales electriques (S U ...,S U S n ) qui sont repartics sur tout le temtoire de Fetal et 

apparaitront figures 3 a 8. Le recepteur (Bo) de mesures observables est sch6matis£ par une 
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antenne parabolique de reception de donnees obscrvablcs (8) qui regoit en provenance du 
satellite (6), en (lux continu pdriodique, ^les r mesures d'environnement 
(M(f t )) = (/i A (fc) ft " L " r) ) dmises par le fournisseur de donnees (Fo) via le reseau (Ro) 
d'emission de donnees observables. L'interface dc modelisation (Ct) est prdftrcntiellcment 
5 silu6e sur Ic site centralist emetteur (U). Dans le cas plus g£n6rai, eile est situee au sein d'un 
site expert (£/')> qui pcut 6tre distinct du site emetteur (U) mais en tout cas distant (de 
plusieurs kilometres, voire de plusieurs centaincs de kilometres) d'unc majority des n 
systfcmes industriels locaux distants (Sj um } '" n) ). Au niveau de ('interface de modelisation (ft) 
ou du site expert (£/')» on fait le choix : 
10- des p variables (X* = 1 " p) ) de renvironnemenl (£), Sgalement appel^es facteurs de risque, 

- des horizons de temps (t k (k s 1 - m) ) des differentes simulations, 

- du module numdrique descriptif d'environnement (d>), definissant ('ensemble de p 

variables d'&al (XAtf mXm '^) - (0>(/i,(/) fi r (t)j) en fonction des mesures observables 

(Mh(0 iho} ' r) ) (avec g£n6ralcmcnt p < r mais suffisamment grand, en sorlc que Pcnsemblc 

15 des variables d'environnement (Xr(/) (l " , "* ) ) decrivent physiquement, & chaquc instant 
d'cchantillonnage (/*), la globalit£ de Petal dc renvironnemenl (E)), 

- et d'une loi statistique param6trte (modele) (fcPC u ...J( p \ /; A/*)) du comportement joint 
des facteurs de risque, au temps (/) el par tranches dc temps (A/** s L "'" ) ) > des variables 
ductal (Xi(0° m} "*\ decrite par sa densitd de distribution jointe. 

20 La loi statistique (f„) param<Slr£e (decrite num^riquemcnt) est gdndralcmcnt chotsie par un 
mathematicien spdcialiste des modeles statistiqucs. Elle estjustifiee 5 partir d'obscrvalions 
historiques — il pcut s'agir notamment de lois normalcs, log-normales, etc., ou plus 
£dn6ralement d-unc distribution statistique jointe des p variables d'environnement (Xfa) 0 = 
- p) ) ou de leurs variations temporellcs (AX/00 °' " "+\ validee par des tests statistiques 

25 classiques adequate (tnoyenncs, variances, correlations, sauls, comportements extremes, etc.). 

Au fur ct a mcsurc dc la reception des donnees observables d'environnement que Ton regoit 
sur demande ou en continu, on met a jour et on enrcgistre dans la mdmoirc dc masse 
(//) Thistorique periodiquc des valeurs des mesures (/i/,(/) (fta ] ~ r) )„ Lc calculateur (C) calculc, 
a chaque periode d'cchantillonnage (/*) successive, h Paidc du modele num£rique descriptif 

30 d'environnement (<[)) et du vecteur de mesures (M(/)), les p variables d'etat 
(XAi)"'* " / ' > ) = (<I>(/ii(/) > ...^r(0))- Prefereruiellement, lc calculateur (O calculc ct met 
egalement h jour nSpetitivemenl Thistorique des variations des p variables d'etat 
(AXfO) 0 ' 1 ~ p) ) qu'il enrcgistre dans la mcmoirc dc masse (#). On entend par variations des 
variables d'etat, soit la difference, soit le ratio de deux valeurs d'unc variable h des instants 

35 successifs, ou plus gdndralement toute valeur numerique d6crivanl les variations (absolues ou 
relatives) de la variable entre deux instants successifs. 

Le calculateur (O est muni d'un systfcme simulateur (9) de scenarios, d'un systeme calibreur 
(10) dc la loi statistique parametree (f„) pour quelle soit conformc aux observations 
historiques memorisces dans la memoirc (//), ct d'un gdndratcur aldaloire (II) permettant dc 
40construire nSpetitivement des « smp series dc simulations probabilistes » (W<f) muhi-valuces 
de description des pseudo-6tats ou pseudo- variations probables futur(e)s des variables d'etat 
(*,</)" s| -~ p) ). 

Diverses melhodes dc fonctionnement du simulateur (9) dc scdnarios seront decrites plus loin. 
Ccrtaincs generent des series de simulations des variables d'etat '(Xftf' l " p \ et d'autrcs des 
45 s6ries de simulation des variations des p variables cTAal (AXfOf* 1 "*'°). En sorte qu'au niveau 
du site centralise dmetteur (CO, a partir de la densitc dc distribution jointe de la loi statistique 
parametree (f a (X u ,,. t X p \ t\ A/*)) retentie du comportement joint dans le temps (r) des 
variables d^tat (X/), on gSnfcre numeriquement repetitivement a des instants / 0 successifs, une 
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« snip sdrie dc simulation probabiliste* (W// <'-'-' , ' sl ■ ■***■»■•»>) multi-valiidc dc description 
des pseudo-etats ou pseudo-variations probables futur(c)s de I'cnvironnement (£>. 
Selon un premier mode de mise en oeuvre, cette «smp sdrie de simulation probabiliste » est 
une « smp sdrie probabiliste » (X t f) de pseudo-etats formde dc sxm /duplets d'dchantillons 
5 muiti-valuds de pseudo-etats futurs (X^ihf^"^^ 1 "^^ 1 *""°), comprcnant sxmxp yalcurs 
numdriqucs possibles des pscudo-valeurs futures probables des p variables d'etat (XX0 (im 1 *' 
et constitutes : 

- pourchacun des m (w > 1) horizons de temps , - w) ) de la simulation, 

. - dc * (avee 5 ires grand » 1) pseudo-etats futurs possibles (X/ V=, "" f) ) = 
10 (XiJtk) iM "+ Jml '" g *) du vecteur des /; variables d'etat (X,(/*) (/ s '■■ ' J) ) . 

Selon un second mode de mise en ceuvre ddcrit figure 5, cette « snip sdrie probabiliste 
variationnelle » (AX//) est formde de sxm /vuplcts d'cchantillons multi-values dc pseudo- 
variations (A^//*y foT " ,|, ^" l " J,l!ai, '" w> ) futures & partir de.la valeur courante des variables ddtat 
(Xi(t o ) 0st , "'' ) ) 1 comprcnant sxmxp valeurs numdriques de pseudo-variations futures probables 
1 5 des /; variables d'etat (X^f* l ~ p) ) et constitutes : 

- pour chacun des m (m > 1) horizons de temps (/*'* B 1 " m) ) de la simulation, 

- de s (avee .9 trds grand » 1) pseudo-variations numdriques futures possibles (AX,-,* 0 = 1 *°) 
= {AX s Jt k ) tM ''+ Jml " J >) du vecteur des p variables d'etat " i " p) ) • 

Comme cela apparait en reference aux figures 2, 3, 4 et 5, grace h I'dmcttcur (Ap) de sdrie dc 
20 simulations probabilistes equipant le site centralist dmcltcur (L0, repetitivement, a des 
instants / u successifs, on diffuse ladite «smp sdrie probabiliste >> (W,/) multi-valude de 
description des pseudo-etats ou pseudo-variations probables futur(e)s, comprenant les sxmxp 
pseudo-valeurs numdriques constitute par : (cas 1) ladhe « snip sdrie probabiliste » (Xif) des 
pseudo-valcurs des variables d'etat (Xijihf^''^^ 1 ' ^^ ' ^), ci/ou (cas 2) ladite « snip sdrie 
25 probabiliste variationnelle » (AX,./) des pseudo-variations (AX u (( k ) {i=} ' " p,J=u " s - hz} ' jn) ). 

Comme on le constate en reference aux figures 2, 3, 4 et 5, cette emission dc simulations 
probabilistes a lieu de maniere ddccntralisee du site dmcttcur (If) en direction de Fenscmble 
dc la multitude des systemes locaux distants (S u ) via le reseau (Rp) d'emission de donnees 
probabilistes. Dans le cas represenie, Temetteur (Ap) de scries dc simulations probabilistes est 

30 schematise par une antenne parabolique (12) d'dmission de donndes probabilistes transmcllant 
successivement les donnees probabilistes, h savoir les « snip scries probabilistes » (W,/) an 
satellite (6). C'cst-&-dire que les rdseaux (Ro) et (Rp) sont confondus. On comprendra, que 
Ton pourrait tout aussi bien utiliser un autre satellite pour remission dc donndes probabilistes 
que eclui utilise par la reception de donndes observablcs, ou un autre type de transmission de 

35 denudes, mitammcnl via le reseau Internet. 

En reference & la figure 3, on voit Torganisation gendrale, selon Tinvcntion, d'un ensemble de 
systemes industries locaux (S u ). II s'agit, dans Tcxemple retcnu lie h Pindustric electrique, 

d'un ensemble de centrales eiectriques (S ]t ...,S u S n ) dont le comportement industricl 

(production dc KWs cnergetiques) est influence par Tenvironncmcnt (E) ddcrit 2l la figure 1. 

40Chaque systeme industriel local (S u iu=] ' M> ) comprend un systeme informatique local (13, 14, 
15) sensiblcmcnt equivalent. L'architecture d'un systtmc informatique local (15) du systeme 
industriel (S u ) est decrite figure 3. Celui-ci comprend principalcment : 

- une interface locale de moderation (Cl u ) pour configurer le site (S„ ( " = 1 n) ) en fonction de 
different* variables et modelcs, 

45 - un caiculateur local (C») pour effectuer les calculs numdriques el statistiques, 

- une mdmoire de masse locale (H u ), 

- un recepteur local (Bp u ) dc donnees probabilistes. concu et rdgle dc manifcre 5 recevoir le 
flux de donnees probabilistes en provenance du satellite (6). 
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On constate que chaque systeme industriel local (S u ) est cquipd d'une antcnne parabolique 
locale (16,17,18) de reception de donnScs probabilistcs cl d'un rgcepteur local de donndes 
probabilistes (Bp u ). De mani&re independante, au scin de chacun des systemes industriels 
locaux (S u ), r6p6titivement on collecte, grSce a I'aniennc parabolique (18) et au r^cepteur 
5 local (Bp u ) y la succession des «smp series probabilistes » (IV//) multi-valudes de description 
des pseudo-dtats ou pseudo-variations probables futur(e)s & sxmxf> valeurs numeriques 
(Wif = X,/ ou AX//) transmises par le satellite (6). 

Par I'intermediaire dc interface locale (Q u ), on choisit w/ (m'< m) horizons de temps, 
retenus parmi les m dils horizons dc temps .(//*" , " m) ). Eventuellement (dans le cas (2) decrit 
10 ci-dessus et repr6sent6 figure 5) ou 1c site industriel local (S tt ) recoit du satellite (6) unc « smp 
serie probabiliste variationnelle » (AX,/), on calculc localement ladite «smp s£ric probabiliste 
des pseudo-6tats » (X//) constitute de pseudo-dchanlillons d'etats fulurs 
(Xi/ti i ) <tml '" PtJm} "'* M * a} " M) ) correspondante (au moins pour les m' horizons dc temps (/*) 
retenus). 

15 Pour chacun des m' horizons de temps (/*) retenus, on applique numcriqucmcnt le modele 
systemique local (D u ) prealablement fixe afin dc d6terminer dans chacun des s scenarios de 
simulation re^us 0' = I ... .*) Tindicateur dc rendement (/„/ = D U (X U1 (X,-//*)) : 

- pour au moins une valeur particuliere (K u ) = (A^i, A„ f/ ) du vectcur (X u ) des paramfctres 

locaux d'action (que Ton pourra etre amene a faire varicr), 
20- et pour chacun des s p-uplets simutes de variables d'etat (X^W 881 ' 7>) ), corrcspondant 
audit horizon de temps (/*) retenu de la serie probabiliste de pseudo-6lat (X,-,//*) 0 a L ,0 ). 

On determine ainsi s pseudo- valeurs simutees corrcspondantes de Tindicalcur de rendement 
(LJh\i)) = D u (\uXXij{h) (,7l] ' p ' J::L[ ''' s) ). On regroupe et on ordonne les 5 valeurs 
correspondanles de Tindicateur de rendement (I u j(t k Xuf s 1 "°)> scion le niveau dc rendement 

25 On calcule et on memorise dans la memoire dc masse locale (H u ) I'histogrammc (F(/ w )) 
ainsi que la distribution dc probability (en nombre ou pourccnlage) (¥„(/„)) du niveau atteint 
de Pindicatcur dc rendement local (/„) dans la configuration (X M ), scion lc modele local (D M ). 
On cn deduit, par un calcul statistique, la valcur d'au moins un paramctre slatistiquc de 
reglagc choisi (V(k u )) Ocl que: moyenne, e*cart type, Value at Risk, utilite, etc..) 

30 caracierisiique de la distribution de probability QV(\ U J)) des nivcaux de rendement (l u 0kX u ) % 

en fonction du vecieur des paramctrcs locaux diction (X u ) = (A„i A w</ ). 

Eventuellement, on effectue localement, au scin de chaque systcme local (5"„), une 
optimisation numcrique (scion des criteres predefinis) du (ou des) parametre(s) statistique(s) 
de rdglage choisi (V(X U )). Puis on r£gle le niveau du vectcur des parametres locaux d'action 

35 (X u ) 5s (A W |, ...,A,, f/ ) de chaque systcme local (S u ) a une valeur optimalc calcule*c en fonction 
des criteres prcctablis, afin notammcnt d'anticipcr Jc comportemcnt probable futur et 
!'efficacit< de chaque systcme (5„) soumis a Tenvironnement (£) au cours du ou des 
horizon(s) de temps (/*) retenu(s). 

Un mode de mise en ocuvre particuliere de T invention decrit figures 3 et 5 vise a rendre 
40adaptaiif pour chaque systeme industriel local (5 J/ (wss| n) ) la moderation de Tiridicateur de 
rendement correspondant a chaque parameire local d'aclion (A tt/ ,). A cet effet, au niveau du 
site expert ((/')> on associe a chaque paramctre local d'action U,,,/'' * 1 * ^ un code generique 
(U> (ne dependant pas du syslemc local (5J). On ctablil pour chaque code (L h ) y et on 
actualise regulieremcnt, un «dvaluateur de rendement » (n fl ) y e'est-a-dire un programme 
45 (materialise sous forme d'un fichicr numeYise) pcrmcttanl de calculer, en fonction du />-uplet 
des variables d'dtat (X) = (X/' = [ " p) ) ct du paramctre local d'action (X uh l la composante 
correspondante du rendement (I uh ) = (n*(^,X)), de telle sorte que, pour chaque systeme 
local (S u ), ledit modele systemique local (DJ consisle simplcmcnt k addilionncr les 
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composantes (/„„) de I'indicateur de rendement correspondant a chaque paramfctre local 
d' act ion (A B /,) : 

^a„.x) = Jn»(^„x) (*) 

Le site expert (IT) est relid au site centralist emctteur par un rescau (R*) et lui transmet les 
5 definitions numeriques adaplatives des composantes de rendement (/„,,/) = (nW/U.X.*)). 
Au niveau du site centralise" emetteur (£/), d'une part, on gdnere numcnqucmcnt et on emcl, 
rdpctitivement a des instants (t 0 ) successifs, a partir du mod&ie statistique rctcnu 
(fa(Xu---J(p, r, A/))du comportemcnt joint dans lc temps (r) des variables d'&at (A',- ( ' = l " p \ la 
«smp serie probabiliste » (Hfy*) multi-valuee de description des pseudo-dtals ou pseudo- 
10 variations probables futur(e)s. D'autre part, a 1' apparition de chaque nouveau code (L h ) ou de 
chaque modification de I'evaluateur de rendement (II/,), on emet de maniere centralists la 
pairc (L h ,U h ) sous forme d'un fichier numcYise\ 

Au niveau de chaque syslfeme industriel local distant (5„) on calcule la composante dc 
rendement (/„„/) = (n„(-l„/„X/)) : 

15 - pour chacun des /»'(/«'< m ) horizons de temps (/<) de la simulation retenus a Techelon 
local (S„), 

" P °"o-5 1 , ia o un des -/'= 1 ( ave< = s tres grand » 1) pscudo-dtats numeriques possibles 

(X, -* ) rccus, ou reconstituts k partir de la serie des pseudo-variations (AX.f), 
- et pour chacun des h = I q paramfetres locaux d'action (X uh ). 

20 Puis, par sommation, selon la formule ci-dessus (*) de I'indicateur dc rendement (/„/), on 
calcule son histogramme, puis le paramfclre statistique retenu (V(\„)). Enfin on applique le 
vectcurs de parametres locaux d'action (K u ) qui optimise la valeur du parametrc statistique 

Un cas particulier du mode de mise en ccuvre ci-dessus consistc en ce que les cvaluateurs dc 
25 rendement (ll,,) soicnl proportionnels au paramfetre local d'action (X uh ) e'est-a-dire que : 

n/,(A„»,x) = x„i, tt/,(X) 

ou (n h ) est appclc « cvaluateur de rendement marginal ». Dans ce cas. le modfclc systemiquc 
local (£>„) revienl a faire un produit scalaire entrc le vectcur de parametres locaux d'action 
(Ku) cl le vectcur de rendements marginaux (n(X)) = (n h (X) ,h ° 1 "•**) : 

30 D„(X U ,X) = (X u | n(X)) (**) 

Lc calcul, au niveau du site local (.?„), du parametrc statistique (V(X„)) s'en trouve simplific et 
son optimisation acc<£leree par l'enregistremcnt prcalablc, dans la mcmoire dc masse (H u ) du 
calculateur local (C„) des series dc rendements marginaux (n tl (Xj ! ) <h = i = '-•"'>). 

Une premiere variante particulierc de mise en oeuvre de I'invenlion est decrite egalemenl en 
35 reference a la figure 4. Elle conceme lecas oil, au niveau d'un site centralist Emetteur (U), on 
gen^rc numdriqucment rtp&itivement a des instants t 0 successifs, a partir du modcle 
statistique retenu (f^X u . . .JC p - i; Ac)) du comportemcnt joint dans le temps (/) des variables 
d'etat (X,), unc serie probabiliste variationneile (AX,/) formee de sxni p-uplcts d'echantillons 
multi-values dc pseudo-variations (AX„W'='" •••*•*=<•••<»>) d - 6lats fu{urs b parljr de V€m 
40courant des variables d'6tal (X i (to?= , -" ) ). Cettc stric de pseudo-variations comprend x x,,,xp 
valcurs_ numeriques de pseudo-variations futures probables des p variables d'ixax 
(X,(h) ■ p ). Elle est consliiuce pour chacun des m (ni > 1) horizons dc temps (^<* = ' ••'»>) dc 
la simulation, dc s (avec s tres grand » 1) pseudo-variations numeriques possibles 
(AXijiit) ~ f ) des p variables d'dtat (X,(( k f m , "' ,) ). On diffuse de maniere repetitive cctte 
45 «smp sene probabiliste variationneile » (AX,/) des pseudo-variations d'etats (AX,//<)" =I - •* 
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, ...i, - ...^ ce de man j ferc ddcentralisde du site dmctteur (LO en direction de Tensemble de la 
multitude des syst&mcs iocaux distants (SJ. Ceci est cffectud par I'dmettcur (Ap) de series de 
simulations probabilistes grSce h Tantenne parabolique (12) d'emission dc donnees 
probabilities ct via le satellite (6). De manifcre inddpendantc, au niveau de chacun des 
5 systemes locaux (S u )> le recepteur local (Bp u ) de donndes probabilistes collecie de manidre 
repetitive via Tantenne parabolique locale (16,I7J8) de reception dc donnees probabilistes, & 
chaque instant (/ 0 ), et selon une frequence dctermindc, voire & chaque modification dc Tunc 
des variables (X/), la valeur courante des variables d'etat (Xi(l 0 f s] "' p \ B collecte egalement 
ladite «smp sdrie probabiliste variationnelle » (AX,/) des variations d'dtats (AXi,j{Ar k f sX - p ' 
]q/»i..j. ku\..M}y en f onct i on de rincrement de temps At k , dmise par le site central (U). Ces 
donnees sont mdmorisdes h Tinterieur de la memoire dc masse locale (//„). Le calculateur 
local (C„) lvalue localcmcnt ladite « smp sdrie probabiliste » (AT,/) dc pscudo-etais 
(AT^W* 5 ' * ssl -' w >) aux instants futurs (/*) = (r 0 + A/a), corrcspondant & la sdrie 
variationnelle par calcul numdrique incrdmental du type : 

On peut utiliser toute autre formule (addition, exponcnticllc, etc..) faisanl intervenir ladite 
« smp seric probabiliste variationnelle » (AX,/) dmise par le site centralise cmettcur (U) el la 
valeur de Pdtat courant des variables d'diat (X,(tf sl ' p) ). 

Unc methode particulidre de misc en oeuvre de cctle premifere variante est celle decrile plus 
20 haul en reference & la figure 4. Au niveau du site centralise dmetteur (U) y on gdndre 
numcriqucmcnt rcpdtUivcmcnt h des instants / c succcssifs, unc « smp seric probabiliste 
variationnelle » (AX,/) telle que ddcrile ci-dessus. On collecte de manure rdpdlitive et on 
memorise la valeur courante des variables d'etat (XjV 0 f = ] ~ p) ). On diffuse de maniere 
rdpdlitive et ddcentralisde du site emetteur (U) en direction dc rcnscmble dc la multitude des 
25 systemes locaux distants (5 W )> 5 la fois : 

- la «smp sdrie probabiliste variationnelle » (AX,/) des pseudo-variations detal 
(AAr/^) <fe, "^ , -'-*" l -- ) ), 

- et la valeur courante des variables d'etat (X,(/ 0 )'' = '■ 7>> ). 

Unc seconde variante de combinaison particulidrc dc mise en ceuvre de Tinvention est decrite 

30 en reference aux figures 3 et 4. Comme cela apparait figure 1, prcalablcment, et en differents 
points de renvironnement (E), on installe une sdrie de r capteurs 1 - r) ), observateurs de 

Fetal de Penvironnement (£), emettant a chaque instant (/) un ensemble de r mcsures de Petal 
do renvironncmeni (H/,0) (h ~ l '" rl ). Dans le cas decrii figures 1 el 4, le site Emetteur (60 s'est 
abonnd h la fourniture dc donndes d'environncmcnt fournies par un fournisscur (Fo) de 

35 donnees observables qui lui-meme opdre le rescau dc capteurs 1 r) ). Au niveau du site 

conccptcur (£/'), irds distant d'une majoritd des n systemes industriels locaux distants 
(S u " ), on definit prealablement un modele numdrique descriptif d'environnemcnl (<I>), 
definissant r ensemble de p variables d'dtat (X,</)" * L ^) = (<1>(/i »(/),..., j/,<0)) en fonction des 
mcsures (jj h (t) ( = r) ) (avec gdndralement < r), en sorte quel'ensemblc des variables 

40d , cnvironncmcn((X,•(/) ( ' =^, p) ) decrivent physiquement, h chaque instant d'dchamillonnage 
(/a), la globalile de Fetal dc Tenvironnement (£). Mais on comprendra que Ton peul 
egalement effectuer cettc operation directement au niveau du site emetteur (U) y ou pour le 
compte de eclui-ci, notammcnt lorsque le site concepteur (Lf) est situd au scin du site 
cmettcur (II). Au niveau du site centralisd emetteur (If), on collecte & distance et on mdmorise 

45 1es r mesures d'environnement Qj h (t) ihss x ~ r \ ces /- mesurcs constituant h chaque instant 
d'echanhllonnage (i k ) les composantes d'un vectcur (emporcl de mesure d'environnemcnt 
(M(t k )) = (JJhOk) <h= 1 °). puis on calcule, h chaque periode d'echantillonnage (f k ) successive, h 
Taidc du modele numdrique descriptif d'environncment (<t>) ct du vecteur de mesures (M(/)), 
les p variables d*dtat (X,(/) </!a 1 •») = («)(//,(/) y^t))). 
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L'invenlion recomrnandc deux methodcs particulidres de simulaiion stochastique rdpdtitive 
hidrarchisec par dchantillonnagc de scenarios probables, pour la prediction probabilistc du 
comportement futur ct/ou Poptimisation des paramfctrcs de rendemeni d'une mukilude de 
systemes industriels locaux distants (£,,). 

5 La premiere mdlhode propose par Pinvention de simulation stochastique rdpdtitive 
hidrarchisde par dchantillonnage de scenarios probables consists en combinaison, ft utiliser la 
methode connue dice de Monte-Carlo. Selon cede premidre mcthode, qui est decrile en 
reference ft la figure 2, au niveau du site centralise emettcur (If) on maintient el on met ft jour 
un hislorique des valeurs des p variables d'dtat (Xftf a 1 p) ) dc Penvironnement global 
lOcommun (F) au cours du temps. On gdndrc numdriqucmcnt rdpdtilivemcnt ft des instants r 0 
succcssifs, ft partir : 

- du module statistique retenu (fotX u ...X v \ /; A/)) du comportement joint dans le temps (r) 
des variables d'etat (X,), 

- et de cet historique mdmorisd des valeurs des variables d'etat {XAff ' m 1 " 7>) ), 

15 la «smp sdrie probabiliste » (IV,/) multi-valudc dc description des pseudo-dtats ou pseudo- 
variations probables futur(e)s de Penvironnement (£), selon la mcthode dite de « Monte- 
Carlo ». C'est-a-dirc que Ton procede par tirages pseudo-aleatoires des pseudo-valeurs (Wij) 
grace 5 un gdndratcur aleatoirc (11) de nombrcs situc sur 1c site centralist emettcur (£/), voire 
sur le site concepleur ((/')• On peut utiliser un processus : 

20 - soil « pscudo-aldatoirc » de nature ergodique, 

- soit « ft discrepance faiblc » (suites de Sobol, Hamerslcy, etc.), 

- soit toute autre mcthode produisant des nombres asymptotiquemenl rcparlis scion une 
distribution spdeifide (par exemple uniforme sur On intervalle). 

Et on diffuse ladile « smp sdrie probabiliste » (W,/) multi-valudc dc description des pscudo- 
25 etats ou pseudo-variations probables futur(e)s dc maniere dcccnlralisde du site emettcur (U) 
en direction de Pcnsemblc dc la multitude des systemes locaux distants <S„). Cette mdthode 
de misc en oeuvre est celle qui est recommandee par Pinvention. 

La seconde mcthode proposce par Pinvention dc simulation stochastique rdpdtitive 
hidrarchisde par dchantillonnage de scenarios probables consiste, en combinaison, ft utiliser la 

30 mdthode connue d'analyse historique. C'est unc simplification de la methode Monlc-Carlo 
decrite ci-dessus. Scion ccttc mdthode d'analyse historique, on dmet la multitude 
d'echantillons des variations probables futures des facleurs de risque directement ft partir de 
Thisioriquc mdmorisd, Au niveau du site centralise dmetteur (U) y on maintient et on met ft 
jour un historique des valeurs des p variables d'etat (W = ] " p) ) dc Penvironnement global 

35commun (£) au cours du temps. On gendre numdriqucment rdpdtitivement ft des instants /„ 
succcssifs, ft partir : 

- du modcle statislique retenu (/•«(*,,... /; At)) du comportement joint dans le temps (/) 
des variables d'etat (AT,), 

- et dc cet historique mdmorisd des valeurs des p variables d'elat (Xj(t) (isi u -+\ 

40 1a «smp serie probabilistc » (W { f) multi-valude de description des pseudo-dials ou pseudo- 
variations probables futur(e)s de Penvironnement (£). Selon la methode dite « historique », on 
reconstruit Pestimatton d' etats futurs probables des p variables d'etat (X&jf= y "**) ft partir de 
Petal courant (X;0 o f= L p) ) et de variations historiqucs passdes (AX,(rO ( ' = -/>) ) pour un 
ensemble de dales passdes (/'< r 0 ). et en affeclant a chaquc etat reconstruit unc probabilitd de 

45 realisation pouvant dependre dc PSge de la variation historique correspondante. On diffuse 
ladite «smp serie probabiliste » (Wj) multi-valudc dc description des pseudo-dtats ou 
pseudo-variations probables futur(e)s, de maniere ddcentralisde du site dmetteur (U) en 
direction dc r ensemble de la multitude des systdmes locaux distants (S w ). 

Une troisidme varianle de mise en oeuvre de Pinvention consiste ft introduire des biais dans les 
50 simulations stochastiqucs cmiscs. Ces biais visent ft prendre en compte des situations 
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extremes qui n'arrivent que rarement et/ou constitueni des « scenarios catastrophes ». Selon 
ceite methode, on biaise la gdndration des «smp scries probabilistcs » (IV//) mulli-valuees de 
description des pseudo-etats ou pseudo-variations probables fulur(e)s dmises par le site 
dmetteur (60 en direction des multiples systdmes industricls locaux distanls (S u ), ce par 
5 rapport k la stalistiquc de Thistorique des variables d*dtats ou des variations dc variables 
d'etat. Selon cette variantc de Tinvention, au niveau du site centralise dmetteur (60* on 
maintient et on met a jour un historique des valcurs des p variables d'dtat (X,</) </= L p) ) de 
Tenvironnement global commun (E) au cours du temps. On gdncre numdriqucmenl 
rdpdtitivement a des instants / 0 successifs la «smp sdric probabiliste » (W,/) multi-valuee de 
10 description des pscudo-dtats ou pseudo-variations probables fu(ur(c)s de Tenvironnement (E), 
5 partir : 

- du modele statislique retenu (f 0 pC\ % ...JC p \ /; A/)) du comportcmcnt joint dans le temps (f) 
des variables d'etat (X/), 

- et de cet historique mdrnorisd des valeurs des p variables d'etat (XM a * l ~ p) y 

15 Mais on inlroduit volontairement dans la «smp serie probabiliste » (IV//) des biais de 
comportement statistique joint des variables d'etat (Xfa) am 1 ,7,) ) vis-ii-vis dc lcur 
comportcmcnt statistique historique. Lcs series biaisees dmises rcproduisent ainsi des 
« scenarios catastrophes ». Ceci permet dc tcnir compte de Timpact desdits scenarios 
catastrophes sur les systdmes distants (S u ). Et on diffuse Iaditc «smp sdrie 

20 probabiliste » (W//) multi-valuee biaisee de description des pseudo-dtats ou pseudo-variations 
probables futur(c)s de manure decentraliscc du site emetteur (60 en direction de Tensemble 
dc la multitude des systdmes locaux distants (S w ). 

Selon unc premidre sous-variante dc la troisidme variance decrile ci-dessus, Tinvention 
rccommande d'adapter selcctivement la leneur des « snip scries probabilistes » (W,/) multi- 

25 valuccs biaisees de description des pseudo-dtats ou -pseudo- variations probables futur(e)s, 
dmises par le site dmetteur (60 en direction de r ensemble de la multitude des systemes 
locaux distants (£„), ce en fonction de requites cffccludcs sdlectivcmcnt par lcs diffcrcnts 
systemes distanls (S u ). A cet effet, certains des systdmes distants (S tt ) envoient au systeme 
centralist emetteur (U) des rcqudtes sous forme de paramdtres numdriques de scdnarios 

30 catastrophes selon des biais specifiques definis par cux. On gdndre au niveau du sile centralist 
(60, en fonction de ces requetes, eventuellemcnt en fusionnant celles qui sont identiques, des 
« snip series probabilistes » (W*/) mulli-valuccs dc description des pseudo-dtats ou pseudo- 
variations probables futur(e)s, biaisdes en fonction des paramdtres specifics rcquis par les 
systemes distants. Le site emetteur (60 envoic selcctivement vers chaque systeme local distant 

35 (S u ), ladite « smp serie probabiliste » (W</) multi-valuee biaisde de description des pseudo- 
dtats ou pseudo-variations probables futur(e)s qui correspond & sa rcqucte, Cet envoi est 
effectud de manidre elccironique (via Internet, un rescau local, le telephone, ou tout autre 
technique permctiant la transmission de donndes numdrisdes). 

Selon unc dcuxidmc sous-variante de la troisidme variante decrite ci-dessus, r invention 
40 recommande que les systdmes distants (S u ) alterent eux-meme lout ou partie des « smp sdrics 
probabilistes » (IV//) qu'ils resoivent du site emetteur (60 pour prendre en compte des 
situation possibles extremes qui les concernent. Selon cette sous-variante, certains des 
systemes distants (S u ) altdrent eux-mdmes tout ou panic des «smp sdrics probabilistes » (IV,/) 
multi-valuccs dc description des pseudo-dtats ou pseudo-variations probables futur(e)s. Ceci 
45 permet dc construirc localcmcnt des sdries biaisdes reproduisant des scdnarios catastrophes 
dont les paramdtres sont : soit specifids au niveau du systeme emetteur central (ZD et diffuses 
en direction des systemes locaux (S u ), soit definis au niveau des systemes locaux (S u ). 

Une quatridme variante de mise en ceuvre prefdree de ('invention consiste en ce qu'au niveau 
du site centralisd emetteur (60, on gdndre numeriquemcnt la « smp sdrie probabiliste » (IV//) 
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multi-valute de description des pseudo-ttats ou pscudo-variat ions probables futur(e)s de 
Fenvironnement (£), de manitre ptriodique, scion une frequence d'echantillonnagc (qf). 

Scion une cinquidmc variantc dc mise en ocuvrc, 1'invcnlion rccommandc trois mtthodes 
pour effeelucr la transmission des «smp series probabilistes » (W 7 ,/) multi-valutes de 
5 description des pseudo-etats ou pseudo-variations probables futur(c)s, du site centralist 
tmetteur (U) vers la multitude des sites industriels locaux distants (S u iu = u ■ n) ) . 

Une premiere sous-variante d'tmission de donntcs probabilistes est decrite figure 6. Au 
niveau du site centralist tmetteur (L0, on diffuse ladite «smp strie probabilistc » (W,/) multi- 
valute de description des pseudo-ttats ou pseudo-variations probables futur(e)s, de manitrc 

lOdecentralisee du site tmetteur (U) en direction de Pensemble de la multitude des systemes 
locaux distants (S u ) par Emission sous unc forme « broadcast », notamment par tmission 
radio. Dc manitre independante, chacun des systdmcs locaux (S u ), majoritaircmcnt trts 
distants du site centralist tmetteur (L0> collecte rtpttitivement ladite «smp strie 
probabiliste » (W}./) multi-valute de description des pseudo-ttats ou pscudo- variations 

15 probables futur(e)s h sxtrixp valeurs numeriques (W,/ = AT,/ ou AX,,/), par reception 
quasiment en continu dc cctte Emission « broadcast » en provenance du site centralist 
tmetteur (t/). 

Une deuxitme sous-variante d'emission de donntcs probabilistes est decrite figure 7. Au 
niveau du site centralise tmetteur (£/), on diffuse ladite «smp stric probabiliste » (W,/) multi- 

20 valute de description des pseudo-tiats ou pscudo- variations probables futur(e)s, de maniere 
dcccntraliste du site emetteur (U) en direction de Fenscmblc dc la multitude des systfemes 
locaux distants (S„) par tmission sous une forme dite « PP aval », c'est-iVdire par tmission 
sur un rtseau (R) —notamment de type Internet — du site centralise tmetteur (£/) a Padresse 
sptcifique de chacun des systtmes industriels locaux distants (S,/" =I n) ). A cct effet, on 

25 constate que le site emetteur (U) est tquipt d % un modem central (20) et que chaque syst&me 
industriel local distant (S„ ( " = 1 - n) ) est tquipt d'un modem local (21) par 1 esq u els transitent les 
donntcs probabilistes. Dc maniere independante, chacun des systemes locaux (S w ), 
majoritairemeni ties distants du site centralist tmetteur (J7), collecte rtpttitivemem ladite 
«smp strie probabiliste » (W,/) multi-valute dc description des pseudo-ttats ou pseudo- 

30 variations probables futur(e)s b sxmxp valeurs numtriques (IV,/ = Xj ou AX,/), par 
enrcgistrement successifs dans la mtmoire de masse locale (H u ) du calculateur local (C w ). 

Une troisieme sous-variante d'tmission de donntcs probabilistes est dtcrite figure 8. Au 
niveau du site centralist tmetteur (t/)» on diffuse ladite « snip strie probabilistc » (IV//) multi- 
valute dc description des pseudo-etats ou pseudo-variations probables fuiur(e)s, dc maniere 

35 dtcemralis^c du site tmetteur (U) en direction de Pensemble dc la multitude des systemes 
locaux distants par mise & disposition de ladite «smp strie probabilistc »(W,/) multi- 
valute sur un servcur (Z) de donntes probabilistes. Cclui-ci est disponiblc sur un rtseau (/?) 
d'emission (notamment de type Internet), ce par acc&s « PP » de chacun des systemes 
industriels locaux distants tsf" De manitre independante, chacun des systfcmcs locaux 

40 (S u ) collecte rcpctitivcmcnt ladite «smp strie probabilistc » (IV,/) multi-valute dc description 
des pseudo-ttats ou pseudo-variations probables fulur(e)s, par acc&s sous une forme dite « PP 
amont», Cest-a-dire quails accedenl au serveur (Z) du site centralist tmetteur (U) via le 
rtseau (R) par requete in son adresse sptcifique (notamment dc type Internet), Puis ils 
acquierent les «smp stries probabilistes » par transfer! (notamment de type «FTP») via le 

45 rtseau (R) et par enrcgistrement dans la mtmoire dc masse locale (H u ) dc leur calculateur 
local (C„). 

Selon une sixieme variante dc mise en oeuvre, ['invention recommande trois mtthodes pour 
organiser des « smp series probabilistes » (Wj/) transmises par lc site centralist tmetteur (U) 
vers la multitude des systemes industriels locaux distants (S/ w = { - "\ ccci afin de compresser 
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le volume des donnecs probabilistes transmises el d'accclerer leur prise en compte par la 
multitude des systdmes industriels locaux distants (5,/ M = 1 '°). 

Selon une premiere sous-variante <P organisation de remission des donndes probabilistes, au 
niveau du site centralisd dmetteur (60, on gdndre numcriquement repdtitivement & des instants 
5 / 0 successifs, & partir du modele statistique retenu (fc6X l ,...,X f> \ t\ Ar)) du componement joint 
dans le temps (;) des variables d'dtat (X^ x " p \ une «smp serie probabiliste » (W f f) multi- 
valude de description des pseudo-etats ou pseudo-variations probables fu(ur(e)s dc 
Tenvironncment (E). Cellc-ci est constitute de s (avec s trfes grand » 1) pscudo-valeurs 
numdriqucs des /> composantes (W^f* 1 ~* ^) pour chacun des m (m > i) horizons de 

10 temps 0k iksU " m) ) t sous la forme d'un «smp fichier » de sxmxp donndes alpha-numdriqucs 
successives . Mais, d'une part, on subdivise ce «smp fichier » en m « sous-fichiers 
temporcls >> successifs de niveau 1 regroupant les pseudo-valeurs d'un mcmc horizon de 
temps Ok s *"'). D'autre pari, chacun des « sous-fichiers tcmporels » est constitud d'une 
succession de j « sous-fichiers dc tirage » de niveau 2 constituds chacun d'un /duplet dc 

1 5 pscudo-valeurs. 

Selon une dcuxidme sous-variante d'organisation dc remission des donnees probabilistes au 
niveau du site centralist dmetteur (60, on gdndre numdriqucment repdtitivement h des instants 
t 0 successifs une « smp sdrie probabiliste » (W,f) multi-valude dc description des pseudo-etats 
ou pseudo-variations probables futur(e)s de renvironncmeni (E) en sdparant les m « sous- 
20 fichiers temporels » successifs de niveau I regroupant les pseudo-valeurs d'un mdmc horizon 
de lemps (/,.<*" 1 ■"°) par un indicateur d'horizon dc temps (rj k iks!} ' ,n) ) corrcspondant. Et, au 
niveau dc chacun des systdmcs locaux (S„), on distingue les uns des autres les different dits 
sous-fichiers lemporcls d'unc mfime «smp sdrie probabiliste » (W,;*) multi-valude collectde 
en identifiant leurs indicateurs specifiques d'horizon dc temps (rfc f *'"' ) ). 

25 Selon une troisidme sous-variante d'organisation dc remission des donnecs probabilistes au 
niveau du site centralisd dmetteur (60, on gdndre numdriqucment repdtitivement a des instants 
f 0 successifs les « smp sdries probabilistes » (W t f) multi-valudcs de description des pseudo- 
etats ou pseudo-variations probables futur(c)s de renvironnement (£}, en sdparant les s 
« sous-fichiers de lirage » successifs de niveau 2 (d'un mdme fichier temporel) regroupant 

30 chacun un />-uplet de pscudo-valeurs par un indicateur de tirage (t/^ 1 -°) correspondant. Et, 
au niveau de chacun des systdmes locaux (S u ) y on distingue les uns des autres les different 
« sous-fichiers dc tirage » successifs de niveau 2 (d'un mdme fichier temporel) regroupant 
chacun un /?-up!ct de pseudo-valeurs d'une mamc « smp sdrie probabiliste » (Wif) multi- 
valude collectde, cn identifiant Ieur indicateur specifiquc de tirage (jxp r = 1 ■ f) ), 

35 Scion une scplidme variante de mise en oeuvrc, V invent ion rccommande trois methodes pour 
filtrer sclcctivement le contenu des «smp sdries probabilistes » (Wtf) transmises par le site 
centralise dmetteur (60 vers la multitude des systdmes industriels locaux distants (S u iu = , 
Ccci vise h adapter Jeur contenu aux besoins et aux particularitds des systdmes industriels 
locaux distants (S,," 55 " n) ) ainsi qu'aux liens hidrarchiques, industriels ou commerciaux entrc 

40 le site centralise dmetteur (60 et la multitude des systdmes industriels locaux distants 

Une prcmidre sous-variante dc filtration dc remission des donnees probabilistes est decritc 
figure 6. Au niveau du site centralisd dmetteur (60, on diffuse dans sa globalitd ladite « smp 
serie probabiliste » (W,j ) multi-valude de manidrc decentralisde du site dmetteur (60 en 
45 direction de ['ensemble de la multitude des systemcs locaux distants (S u ) par emission sous 
une forme « broadcast ». On peut notamment proedder par emission radio, accessible a tous 
les systemcs Jocaux (S u ). El, de manidrc inddpendantc, au niveau de chacun des systdmes 
locaux (S U X on installe et on utilise un filtre seleclif prereglc (22) de mise 5 disposition d'unc 
fraction (IV,./) preddfinie de ladite «smp scric probabiliste » (Wif) collectde par reception en 
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provenance du site centralise emetteur (IT). Les caracttrisliques de filtrage des donndes par ce 
filtre selectif peuvent notamment ctre acquises par un abonncment spdeifique pris par chaquc 
systtme local (S u ) auprts du site tmetteur (U). Ellcs pcuveni ctre modifiees par ttltaelion h 
travers le rtseau reliant le site centralist emetteur (U) aux systemes locaux (S u ). 

5 Une deuxteme sous-variante de filtration de remission des donntcs probabilistes est dtcrite 
figure 7. Au niveau du site centralise 6metleur (U), on diffuse de manure selective (scion une 
r&gle d'attribulion prddefinie) diverses fractions difftrcntcs (W/)^* 1 • ,) ) prtdefinies 
(notamment acquises par abonncment sptcifique) de ladile «smp scrie probabiliste » (W,,/). 
Cette diffusion selective est effectute de mani&rc dtcenlralisce du site Emetteur (U) en 

10 direction de plusieurs groupes (G M (J}s=l ~ v) ) attribulifs dc la multitude des syst6mcs locaux 
distants (S u ). Chaque groupe (G R ) recoil une fraction particulierc Ccci s'cffcctuc par 

Emission sous une forme dite «PP aval », e'est-a-dire par emission sur un rtseau (R) 
(notamment de type Internet) du site centralise tmcucur (U) h Tadresse spdeifique de chacun 
des systemes industriels locaux distants (5 4r (tfo| -" ) ). De manure independante, chacun des 

15 groupes (G x * m ^ v) ) attributifs dc syst&mes locaux (S u ) collcctc rcpetitivemcnt ladite fraction 
predefmie ((W,/W de la « snip s6rie. probabiliste » (IV,/). lis enregistrcnt successivemcnl 
cette fraction predefinie ((W,/\) dans la mtmoire de masse locale (H u ) de leur calculateur 
local (CJ. 

Une troisteme sous-variantc de filtration de remission des donndes probabilistes est dtcrite 
20 figure 8. Au niveau du site centralist emetteur (U) y on diffuse ladite « snip strie 
probabiliste » (W,/) de manifcre decentralisde du site emetteur (II) en direction de Tensemble 
de la multitude des systemes locaux distants (S u ) par misc ^disposition de ladite « snip strie 
probabiliste »(\V itj ) sur un serveur (2) disponible sur un rcscau (R) (notamment de type 
Internet) auquel chacun des systemes industriels locaux distants (5,/" = 1 *" ) ) petit acceder par 
25 accfcs_<<PPamont>>. Mais on met de manttrc selective h disposition de diffdrcnts groupes 
(Gg** v ) de systemes locaux (S u ) sur le serveur (Z) diverses fractions limit at ives (OV/j \.) 
predtfinies (notamment acquisc par un abonnement sptcifique) dc ladile « snip strie 
probabiliste » (W//). Ceci pcut notamment s'optrer par Pinterrnddiaire d\m sysifcrnc 
d'identificateurs et dc mots dc passe (23) transmis par chaque systemc local (5,/ M5S 1 lors 
30 de son accds au serveur (2) . De manure independante, chacun des systemes locaux (S u ) 
collecte rtpetitivemcnt ladite fraction ({Wif) g ) predefinic attribute dc ladite « smp strie 
probabiliste » (W,-/) par acc£s sous une forme dite « PP aval », e'est-h-dire par acefcs via 1c 
rcscau (R) du site centralist emetteur (U) b son adresse spdeifique. Us acquifcreni par transfert 
ladite fraction ((Wg/) x ) predefinic attribute dc ladite « smp stric probabiliste » (W,/) multi- 
35 valuee qui leur est attribute. lis enregistrcnt successivement cette fraction prtdefinie ((W,/),) 
dans la memoirc dc masse locale (H u ) de leur calculateur local (C). 

Selon une huitidme variante de mise en ceuvrc, dtcrite en rcftrcncc h la figure 10, Pinvention 
recommande d'effectuer une optimisation probabiliste groupte de certains groupes de 
systemes industriels locaux distants (S u (il * 1 % Ceci se rapportc notamment au cas ou les 

40systtmes industriels d'un meme groupe d'optimisation ont des objectifs industriels (de 
production, d'economie, etc.) communs qu'il convicnt d'oplimiser globalcmenl. On choisit 
prtalablement n' [n'<n) systemes industriels locaux, rclenus parmi les n dits systdmes 
industriels locaux (5 W ). On choisit dgalement prtalablement m'(m'<m) horizons de temps, 
retcnus parmi les m dits horizons de temps (//* = 1 ~ M \ identiqucs pour tous les systemes 

45 locaux. 

Au niveau de chacun des systemes locaux (S u ) retcnus, et pour chacun des horizons de 
(emps (u) rctcnus, on applique num£riqucment le modclc systtmiquc local (D„) prtalablement 

fixt afin de determiner dans chacun des s sctnarios (J = 1 .y) lc niveau de I'indicateur de 

rcndcinent (I u f = D u (X uy (X fJ (t k )) : 
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- pour au moins unc valcur particuliere (k u ) du vccieur des parametres iocaux d'action, 

- et pour chacun des s /duplets simules de variables d'etat (X t j(t k ) a=l p \ correspondant 
audit horizon de temps (t k ) retenu, de la serie probabiliste de pseudo-6tals 

5 On determine les s pseudo-valeurs simutees de Tindicateur de rcndcmcnt 
Uuj(tiiXu)) = (DuCkuyXjOi;))^* 1 ■ correspondantcs. On 6mct la « se*rie de simulations locales » 
constitute des s pseudo-valeurs simuldes de Tindicateur de rcndcment 
Uuj{?k*K)) = (/>>«(^«.Xy(/ A .)) 0 = l - correspondantcs du syslcrnc local dc manicre 

central isantc cn direction d'un site centralisateur (B). 

lOAu niveau du site centralisateur (5), on collcctc les sxn' pseudo-valeurs simulees 
d'indicateurs de rendements UujihXu)) 0 * l " J) ) en provenance dc chacun des n' systemes 
locaux (S u ) et pour chacune des s simulations, on calcule, pour chaquc simulation (j = I, ...,5), 
un indicateur de rendement global (JJ) de l'ensemble des systemes locaux & partir d'un modele 
devaluation global (Jj) = (x(!uj(tM iU * 1 ' * n J=1 On regroupe et on ordonne les j valeurs 

15 correspondantes de Findicatcur de rendement (7//* ..,V) 0= l ' ^ selon son niveau. On 
calcule et on memorise l'histogramme (T) ainsi que la distribution dc probability (en nombre 
ou pourcentagc) (H J (f)) du niveau atteint de Tindicatcur dc rcndcmcnt global (J) dans la 
configuration choisiedes vecteurs de parametres locaux d'action (A) = (A,,/" = , -'°), selon les 
modeles locaux (D„) et le modele devaluation global (x) rctcnus. On en deduit, par un calcul 

20 statislique, la valeur d'au moins un paramctre statisliquc dc rcglagc global choisi (V(A)) (tel 
que: moyenne, ecari type, Valuc-at-Risk, utility, etc..) caracteiistique de la distribution de 

probability (x(?h \ n \ J)) des niveaux de rendement globaux (JOk 'Xi *«■)) en fonction 

des valeurs du muiti-vecteur (J) = (Lll,...,Gn() des vecteurs dc parametres locaux d'action 
(Liu). 

25 Eveniucllement on effectue globalcment une optimisation numenque (scion des criteres 
predefinis) du paramctre statistique de rcglagc global choisi (V(.J)) par rapport au multi- 
vecteur (..') = (J l,...,Cln). Et on transmct h chaquc systcme local (Su) 1c niveau 
correspondant de son vecteur optimise des parametres locaux d'action (\ tf ), afin notammcnt 
d'amiciper globalcment 1c comportement probable futur et Tefficacitd de Tcnscmble des 

30systemc locaux (5 W ) soumis a Tenvironnement (E) an cours du (ou des) horizon(s) de temps 
(At) retenu(s). 

On decrit en reference h la figure 9 unc application du procedc scion V invention dans le 
domaine dc Tinduslrie financiere. Le procede de Pinvcntion effectue une simulation 
stochastique repetitive hicYarchisee par dchantillonnagc de pseudo-6lats futurs probables de 

35 variables de marchc, pour la prediction probabiliste du comportement futur et/ou 
r optimisation des parametres dc rendement d'une multitude de n (n etant ires grand » 1) 
d'opcralcurs financiers locaux (S u <uzs l "" ) ). L'inycniion s'applique dans le cas oO les 
operateurs financiers locaux (S u (u = 1 ~ sont decentralises ct trcs distants les uns des aulres. 
Chaque operateur financier local (S u iUB u rt) ) possede un portcfeuille (P u itt=s 1 ' n) ) constitue' 

40d*actifs financiers (d h {h * 1 ••«>), dont la valeur (f u (ttal " n) ) depend des mouvernents du marchd 
des capitaux (E). La valcur (/ H (tt=, "" > ) d'un portcfeuille (Pj^ K ~ n) ) pcul etre calculcc par 
combinaison d'un vecteur de variables de marchc (X(0) = (X,{tf = " p) ) et, pour chaque 
operateur financier local distant (5 P ), d'un vecteur dc quantites (X u ) = (A w/3 ( * = 1 v> ) dc chacun 
des actifs financiers {9 h ih = x detenus dans son portefeuille (P tt ). Les variables de 

45 marchc (X { (t) u = ~' p )soM communes aux differcnts operateurs financiers locaux (S u (u " 1 ■" ) ). 
Elles varicnt au cours du temps (/). n est possible pour chaque operateur financier local 
distant (S u ) de reconstiluer avec une precision acceptable la valcur dc tous les actifs financiers 
pouvant ctrc presents dans son portefeuille (P u ), selon un modele local devaluation (D u ), a 
partir du vecteur des variables de marche" (X(r)) = (Xftf m l " p) ). Les variables dc marcl^ 
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peuvent reprdscnier des prix d'aetifs ou de commodity dcs laux dMmSrftt ou des coefficients 
de courbe dc Laux, des volatility ou des coefficients de surfaces de volatility dcs indicateurs 
de defaut d'un dmetteur ou d'une contrepartie, ou encore tout autre nombrc apportant une 
information pertinente sur Fetal des marches de capitaux. 

5 Par exemple la valcur (/,<*-'■ •">) d'un portefeuille (P u (uts ] - n) ) pcut etre calculi selon la 
combinaison lindaire : 

UO^A^TTj/.XW.e,,) (***) 

00 *>,(/,X(/),0 A ) represent lc prix de l'actif(0,,) a l'instant (/) selon le modele local (D„) 
lorsque le marche est dans Fetal (X(/)), ou une autre formule pouvant faire apparattre une 
lOnon-linearite par rapport aux quantity (A,,,,), pour tenir comptc notammcnt du risque dc 
liquidity. 

Selon une premiere sous-variante duplication du procede scion I'invention au domaine de 
Industrie financiere, au niveau d'un site expert tf/'), on fait prcalablcment lc choix : 
- dcs p variables de marche (*/'" •■•*>) decrivanl P ensemble du marche (£), 

15- e( d'un modele statistique du comportcment joint au temps (/) et par tranches de temps 
(Ah - "■") des variables de marche (X,(tf = i p \ decrit par sa densitc de distribution 

jointe (fo(X X p ; t; At k )). 

Le^modcle statistique du comportement joint au temps (/) ct par tranches de temps (A//* = 
) des variables de marche est justifie $ parlir d'observations historiques — il peut s'agir 

20 notammcnt de lois normales, log-normales, etc., ou plus gendralcmeni d'une distribution 
statistique jointe des p variables d'environnemenl (X,{hf = p) ) ou de leurs variations 
tcmporelles (AX.OO' " ■'' ). II est valide- par des tests statistiques classiques adequats 
(moyenncs, variances, corrdlations, sauts, comportements extremes, etc.). 
Au niveau d'un site centralise cmetteur (U), on genere numcriqucment, repctitivement a dcs 

25 instants (/„) successifs, h partir du modele statistique retenu (f a £Xu...X P \ /; A/ t )) du 
comportement joint dans le temps (/) des variables dc marche (AT," = une «smp serie 
probabilistc » (VV, V ) multi-valuec dc description des pseudo-Slats ou pseudo- variations 
probables futur(e)s du marche (E). 

Scion un premier mode dc mise en ccuvre, cette «smp scrie probabilistc » (W,./) est une « snip 
30 scrie probabilistc » (X,/) de pseudo-etats formee dc sxn, />-uplets d'echantillons multi-values 
de pseudo-etats fulurs (XM' mt - p ' ,mt ~^ tm, "^\ comprenant sx,„xp valcurs numcriqucs 
possibles dcs p variables de marche (X,<0 (,= '-'"). Elle est constituee : 

- pour chacun des m (m > 1 ) horizons de temps = 1 dc la simulation, 

- dc s (avec i tres grand » 1) pseudo-etats future possibles des p variables dc marchd 
35 (X i (h) l '- , " p) ). 

Selon un second mode de mise en ceuvre, cettc « snip serie probabilistc » (H',/)est une « smp 
serie probabiliste variationnelle » (AY,/) formee de sxm /duplets d'echantillons multi-values 
de pscudo-variations_(A^ v {/,) < '= 1 - ' , -> 1 •••"'>) fu(ures a partir de relal courant dcs varjables 
dc marche (X,(t 0 )'- ), comprenant sxmxp valeurs numdriques de pseudo- variations 
40 futures probables des p variables d'dtat (Xfaf - ' ■ Elle est constituee : 

- pour chacun dcs m (m Z I) horizons de temps (r A . (l = '•■•'">) de la simulation, 

- de s (avec j tres grand >> 1) pseudo-variations numcriques possibles (AX,Xr*) v " des 
p variables d'etat (XKttf" 1 "^). 

Lc procede de l'invention consiste alors, au niveau du site centralist dmctteur ((/) tres distant 
45 de la majonte des n operateurs financiers locaux distants (5„ l " = 1 " , ), repetitivement a dcs 
instants 0 o )en ce qu'on diffuse ladite «smp scrie probabiliste >» (VV,/) multi-valucc de 
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description des pseudo-elats ou pseudo variations probables fuiur(e)s comprcnant sxmxp 
pseudo-valeurs num<5riques constitute par: 

- (cas 1) ladite « smp serie probabiliste » (Xif) des pseudo-valeurs des variables de marche 

(X/A) " "O. 

5 - (cas 2) ct/ou ladite « smp sdric probabiliste variationnelle. » (AST//*) des pseudo-variations 

(AXi J [ft) Vml "+< Jmt ■■**-'-■»»).. 
Cettc diffusion de ladite «5»ip serie probabiliste » (W;/) est effectuee de maniere 
decentralisee du site emetleur (60 en direction de 1'ensemblc de la multitude des operateurs 
locaux distants (5„) . 

10 De maniere independante, au niveau de chacun des operateurs locaux (SJ, majoritairement 
tres distants du site centralise emetteur (If), on choisit m\m'< m) horizons dc. temps retenus 
parmi les w dils horizons de temps (// =I - J " ) ). Evcntuellement, dans le cas (2) ou I'on rccoit 
umquement ladite « smp serie probabiliste variationnelle » (AX,/) d'cchantillons variationnels 
(AXiji( k )\ on calcule localemcnl ladite « smp sdrie probabiliste (X,,*) des pseudo-etats » 
15 constitute de pseudo-echantillons d'etats futurs (X;.fi.t k ) u * x *-'••••'<>) correspondante, (au 
moms pour les m' horizons dc temps (t k ) retenus). Puis, pour chacun des m' horizons dc 
temps (/<) retenus, on calcule numeriquement, scion lc module local d evaluation (D„) 
prealablement fixe dans chacun des s scenarios (/ = l,..,v), la valcur globale (/„/) = (D M (X„, 
(*///*)) du portefeuille (/>„). On effectue ce calcul : pour au moins unc valeur particulicrc du 
20 vectcur (X*) = (A,,,,...,^) des quantitts d'actifs detenus ct pour chacun des s p-uplcts simulds 
de variables d'euu (AT///*)), correspondant audit horizon de temps (t k ) rctcnu dc la serie 
probabiliste de pseudo-etat CXtfr*) 0 ' ' " s) ). El on determine les s pseudo-valeurs simulces du 
portefeuille <f„Jt k ,X u )) = (D u a u ,X/t k )^"- •■*>) correspondante.s. On regroupe el on ordonne ces 
j valcurs du portefeuille (/»„). Puis on calcule et on memorise riiistogramme (r„) ainsi que la 
25 distribution de probability (en nombre ou pourcentage) (¥.</,,)) de la valeur (/*) du 
portefeuille local (P u ) dans la configuration (\„), scion lc modele local (£>„). On en deduit. par 
un calcul statistique, la valeur d'au moins un paramfctre statistique dc reglage choisi (V(X U )) 
(tel que : moyennc, ecart type, Value at Risk, utilite, etc...) caracteristique de la distribution 
dc probability QVfajy) des valcurs (l u (i t ,(K u )) cn fonction des variations du vecteur 

30 (k u ) = X ut/ ) decrivant le ponefeuillc (P u ). 

Eventucllcment, on effeclue localement, au niveau dc chaquc operatcur local (S„) une 
opt.m.sat.on numcriquc (selon des critercs preddfinis) du (ou des) paramctrc(s) statistiquc(s) 
dc reglage choisi (V(XJ), et on ajuste la composition du portefeuille (\„) = (A,„, ...,A U< ,) de 
chaque opcrateur local (S„) h une valeur optimale calculee en fonction des critercs preetabiis 
35afin notammcnt d'anticiper le comportement probable futur et le rendement de chaque 
portefeuille (/>„) soumis au marche (E) au cours du ou des horizon(s) de temps (/,) rctenu(s). 
Scion unc deuxieme sous-variantc duplication du precede scion l'invcntion au domaine de 
I Industrie financiere, pour evaluer la valeur (/„/) du portefeuille (/>„) d'un opcraleur financier 
OS,,) a la date future (/*) selon le modele (D„) ct dans Vim du marche Simula (X,</*)) 
40 correspondant au sctnarioy, on prend en compte dans le modele (D u ) les ctats dc march* 
s.mulcs, dans ie memc scenario j, aux dates ulterieurcs (,/""*'••""), scion lc principe suivant : 

- on definit pour chaque actif financier ((0„) ( * = '■••*>) une quantite maximale (v„) qu'il est 
possible dc trailer (h 1' achat ou a la vente) dans une jounce, 

- si aucunc des quantity (A„,,) d'actifs (ft,) detenus dans les portcfeuilles (P„) n'excedc la 
45 valcur l.mile (v h ), alors on applique la formule (***) pour evaluer, a la date (n) et dans la 

simulation (J), la valeur (/„/) du portefeuille (/>„) : 
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- si unc des quamitds (A u/l ) d'actifs (ft,) ddtenus dans un des portcfcuilles (P w ) exefede en 
valeur absolue la valeur limite (v h ) y alors on rempiace la quantite U uh ) par (v h ) % affect^e de 
son signe, & la date (/*) et on reporte 1e reliquat {X uh ± v h ) sur la dale suivante (tt*t) en 
utilisant le vecteur de variables de march<5 (Xj(t k + t )) ; 
5 - si cctlc quantite est encore supdricure h la limite (v h ) en valeur absolue, on la fixe ii la 
valeur maximale (v h ) t avee le m6me signe, et on reporte 5 nouveau le reliquat (X uh ± 2\> h ) 
sur la date suivante (t M ) en utilisant le vecteur dc variables de marche (Xj(t M )), el ainsi 
de suite. 

Cede disposition particultere permet de tenir compte du risque de liquidity imposant des 
lOcontraintes de borne supdrieure a la quantity de chaque actif (ft,) que Ton est en droit de 
valoriser selon le module local (DJ a une date donnee (/*), et entramanl un report evcntucl dc 
la position virtuelle rcsuhante dans chaque simulation sur des dates ultcricurcs. 

Applications industrielles de Pinvention 

L'invention trouve une application induslriellc dans tous les domaines (els que l'industrie 
15 manufacturifcrc, dncrgctiquc, agroalimentaire, financterc, dc transport, de sante publiquc ou 
d'exploitation des matures premieres, ... soumis a un environnement global commun et oti il 
convicnl de mesurer et/ou dc controler P impact sur la production des evolutions possibles des 
variables d'environncmcnt globales et locales cn unc multitude de sites industries eloignds. 

A titre d'exemple, on pent avantageusement appliqucr les enscignemcnts de Pinvention 
20 pour Poptimisation de Parrosage de cultures, e( pour le calcul des temps de transport routicr 
ou la gestion d'une flottc automobile. 

Ainsi, une premiere application de Pinvention concernc P agriculture. Pour cette application, 
les systcmcs industriels locaux distants (Sj Mmt '" n> ) sont les exploitations agricoles dc 
diffcrcntcs parcelles eultivdes, Les variables d^tat (XKi) Uul '" pr ) dc Pcnvironnement sont la 

25 temperature dc Tain la piuviometrie, Phygromctric, la sccheresse du sol, etc. cn differents 

points geographiques. Les parametres locaux d'action (X„) = (*„, \ uq ) son( principalement 

les caracteristiqucs des parcelles et des plantations. Le but de Pinvention est de limiter la 
surconsommation d'eau cn declenchant les diffdrenls systemes d'arrosagc automatiques 
uniquement lorsquc la probability que les cultures aient a souffrir de sccheresse dans les 

30 heures ou les jours a venir est sup^rieure a un seuil d'acccptance defini par Pexploitant dc 
chaque parcelle. Ainsi, les systfcmes d'arrosage automatique ne se declcnchcronl pas si la 
probability d'une averse dans les heures a venir est forte. Invcrsement, ii pourra par exemple 
sc declcncher automatiquemcnt un soir si la probability d'une forte s&rheressc le lendemain 
est elevde. La mutualisation entre les exploitants de la structure de generation d'echantillons 

35 do variables mctcorologiques et done des frais lids au site centraiisd ^metteur probabiliste (U) 
rend la mise cn place de ce syslfcrne de ddclcnchcmcnt d'arrosage automatique moins 
couteusc. 

Unc seconde application de Pinvention concerne Poptimisation de trajet routicr. L'objectif est 
d'estimer la loi dc probability du temps d'atteinte de la destination d'un automobiliste ou d'un 

40 transporter compte tcnu de F ensemble des facteurs aleatoircs cxtyrieurs. Pour un conducleur 
qui doit sc rendre cn un point h une hcurc donnee, sculc Poptimisation d'un critere 
probabiliste complexe (probability d'atteindre le point fixy a Fhcure donnye) pcui lui 
pcrmettre de rycllcment optimiscr Fitineraire k choisir. En effct, slatistiquemcnt, il peut 
arriver qu'un tron^on de route soit gynyralcment fluide, mais que Iorsqu'un ralentissement 

45 survient, celui-ci restc encombre longtemps : la vision du temps de trajet moyen (espere) est 
de ce fait trop reductrice. Selon Pinvention, chacun des vyhicules d v un ryseau routier constitue 
alors un sy^steme industriel local distant (S„). Ceux-ci sont soumis aux facteurs aleatoires 
(X,it) ts - p ) qui sont par exemple : la fluidite de la circulation sur les diffyrents tron<;ons 
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routiers, la situation meteorolbgique, la survenue d'unc panne ou cTun accident... Lcs 
paramfctres locaux d'action (X w ) = (A w i, ...,A u<y ) de chaque vehicule sont la situation sur le 
rcseau routier, Titindraire h venir, la vitesse du vehicule, le type de vehicule et de conduite, 
etc. Au niveau des differents systfemes industriels locaux, Tintcrfacc (£i w ) de saisie des 
5 paramfclres (A W |, ...,A tt<y ) peut &tre relive aux differents instruments dc bord du vehicule 
(notamment h un systfcme de navigation par satellite de type GPS ct des captcurs du vehicule) 
et le r£cepteur (Rp u ) de donnees probabilistes peut etre lie h un terminal Internet embarque. 

Une Evolution dc cette sccondc application conccrnc la gestion d'une flotte automobile 
(transport en commun, transport de fret, loucur de voilure, etc.). Si chaque automobile ou 

lOsystdme roulant est dquip<§ d'un systfcme dc localisation (par exemple de type GPS) cl d'un 
emetteur, il est possible de connaitre la repartition de toute la flotte a un instant passe ou 
present, n est aussi possible de connaTtre la situation «la plus probable » (l'espdrance) de la 
position des vehicules & une date future donndc. En centralisant au niveau d'un site 
centralisateur (B) les pscudo-valcurs simulees des indicalcurs de rcndcmcnl (/.//jtAj), 

15 Tinvention pcrmct au gestionnaire du pare d'affincr sa vision sur la probability des 
ev^nements h venir. Ainsi, le gestionnaire du pare peut controlcr automatiqucmcnt la 
probability de la survenue d'evfcnements peu frequents mais lourds de consequences 
(probability que tous les camions d'un transporteur de fret tombent en panne simultandment, 
probability pour que tous les vdhicules d'un loucur dc voiture se retrouvent concentres dans 

20 une region). 

Une troisifcme application parliculi&rement avaruagcuse de ]' invention concerne la maitrise 
des parametrcs d*unii6s industriellcs fonctionnant sur un m6me reseau et soumises h uh 
environnement global alcaloire commun. A litre d'exemple, on peut avantageusement 
appliqucr lcs cnscignemcnls de l'invenlion au domaine de la production d'clcctricite, dans un 
25 environnement ddrdgule, o& le prix de r electricity est fixe par confrontation de Toffrc dc n 
productcurs (S u (uzs l ' ny ) et dc la demande d'un nombre elevd dc consommateurs sur un marchd 
central (le « pool dnergetique »). La rencontre de I'offrc et de la demande peut se faire sur la 
base d'un planificateur central, indiquant h chaque producteur quels doivent 6tre ses 
parametres de fonctionnement. Une telle approchc pose un double problfcmc : 
30- Le planificateur central doit prendre en compte une masse considerable d'informations, 
non sculcmcnt des informations communes (le niveau de la demande, l'£tat du reseau, lcs 
prix d'encrgic primaire), des informations locales interagissant fortement avec 
I'environnement global (le climat), mais ggalemcnt dc multiples informations spdeifiqucs 
a chaque producteur cl independantes dc Tcnvironnemcnt global (Petat dc marchc et lcs 
35 performances de telle turbine, les probtemes de limite de rejet polluant dans telle rivfere, 
etc ...). Scion lcs methodes classiques, cc problcmc d'optimisation devient insoluble dans 
ccs conditions, sauf & faire des hypotheses grossiererncnt simplificalrices sur lcs 
contraintes locales. 

- Chaque producteur local doit s'efforcer d'etre capable dc repondre 5 la requete du 
40 planificateur central. Ignorant par definition 1c comportcment du module d'optimisalion 
central, il se trouve confronte a un aiea de type administratif (quelle requete va adresser le 
planificateur central ?) en sus d'un environnement global aleatoirc. Cet alca administratif 
est d'autant plus difficile & modeliser quMl ne prend en compte qu'imparfailcmcnt lcs 
contraintes operationncllcs locales. Le risque opcrationncl devient alors difficile h gerer. 
45 Dans ccs conditions, Poptimum global selon lcs methodes classiques devient illusoire, car 
il repose sur unc seric ductals opdrationncls locaux loin d'un reel optimum. 

Cette approche, predominate dans 1c domaine de la production d'clectricite, a ete analyst 
par les demandeurs. II ressort de cette etude qu'elle conduil en pratique h des situations dc 
surcapacites industrielles et energetiques significatives. 
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L'akernativc classique pour ddlocaliser le problfeme d 1 optimisation esi d'etablir un marche de 
rencontre d'offre et de demande, avec un prix d'equilibrc. La theorie dconomique montrc que 
la recherche d'optimums locaux correspond, sous certaines hypotheses trfcs gdndrales de 
convcxite, h un optimum global. En mati&re d'dlectricitd, le probleme est techniqucment ardu, 
5 car relectricile est un bien non stockable. H y a done autant dc marches que de pdriodes de 
Tannde. Le mouvemcnt general recent de dereglementation des marchds dc relectricile 
conduit aujourd'hui a plusieurs modules d'organisation de marche, qu'il serait faslidicux de 
ddcrire ici. Mais ils se rapprochent tous plus ou moins du modele ci-dessous. 

Le «pool encrgdtique » est segmentd en (m) sous-marches, correspondani a des horizons 
10 temps de « prdavis » distincts (A;* <fc = , " m) ). A titrc d'cxemple, on distingue un marchd trfcs 
court terme (quclques minutes) destind a maintenir la tension sur le reseau, un marchd court 
teime (1 heurc, « balancing market ») pour compenser les alcas dc Toffrc et de la demande, 
un marche « normal » a 24 heures, etc... Chaque producleur (Sj? u l (usinc de production 
electrique) doit s'engagcr, pour chaque sous-marchd, sur le respect d'unc courbc de cout dc 
15 production (C* w (kO), donnant un cout de production pour chaque niveau de capacite (*•) offcrt 
au rdscau (ce cout pouvant ctrc infini, au-dela d'un niveau maximum de capacite)! Le « pool 
dnergetique » classe sur chacun des sous-marches les offrcs de capacite par cout croissant et 
les confronte avec la demande sur ce sous-marchc : 1c prix (P k {l " = " r> ) du kilowatt-hcurc sur 
1c sous-marchc resuhe de eclte confrontation. Les producicurs donl le cout de production est 
20supdrieur au prix (P k ) nc sont pas appeids par le gcsiionnaire du rdscau. Si un producteur, suite 
a une indisponibilile accidentelle, constate, un temps (A/*) avani la livraison, son incapacild h 
respecter son engagement sur le sous-marchd 1 (avec At k < An ), ii doit racheter les quantitds 
correspondantes sur le sous-marchd (k) (ou tout autre sous-marchd ayant un prdavis plus 
court). Dans ce cadre duplication de Pinvcntion, I'indicateur de rendement du producleur 

25est Jl a i T rge brutc ' Les P aram& tres de commandes (X u ) = (\ ui A™) d'un producleur 

(^-'• ^sontnotamment : 

- Tallocation de sa capacite disponible (Cap u ) cnlrc les m sous-marches et un w-H-icmc 
sous-marchd spdeifique qu'est l'entretien programmd dc la centralc. 

- le choix de Tcncrgie primaire a utilise^ lorsqu'il a cctte faculte (gaz 7 charbon, fuel / 
30 charbon par excmple). 

Le choix des paramctres est soumis a des conlraintes opdrationnelles locales tel que le temps 
ct lc cout minimum de ddmarrage d'unc unite ou de changement de type dc combustible, les 
minima et maxima techniques, les limites de rejets polluants dans Vair et dans Teau, etc. Le 
probleme d'optimisaiion local est complexe, dans la mesurc oO il fait appel a la prdvision sur 
35 la base de variables alcatoires relevant non seulemenl d'un environncment spdeifique a 
('installation (disponibilite des installations de production), mais aussi d'un environncment 
global commun a tous les producteurs, comprcnant : 

- le climat, et notamment la tempdraturc pour une centrale therm iquc (affectant le 
rendement ihcrmodynamique de la centralc), la pluviometrie pour une centralc 

40 hydraulique (affectant debit et rdserves). . . 

- les prix sur les m sous-marches d'electricite, 

- les prix des combustibles utilisds par le producteur, etc. 

Cet environnemcnt global est bien commun a I'cnsemblc des producteurs : les prix sont les 
mcmcs pour tons et les variables climatiques sont non seulcmcnt fortcment corrdldes entrc les 
45differents siles mais interagissent de surcroit avec les prix (on con?oil par exemplc qu'une 
journec irbs froidc prdvue a 24 heures se traduise par une cnvolee du prix de 1 electricite sur le 
marchd a 24 heures). 

Le probleme d'oplimisation global par un planificateur central est rcmplace par n probldmes 
d'optimisation locaux, afin par exemple de maximiscr Pesperance de marge brute et de 
SOgaranlir a x% la marge brute minimale pour supporter la charge dc la dette (typiquement un 
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calcul de Value-at-Risk), ccci impliquant des calculs complexes. Selon Tan anterieur, les 
solutions possibles sont les trois suivanies : 

- Chaquc producteur contracts aprfcs mise en concurrence, avec un organismc (« asset 
manager ») disposant des outils d'analyse de Penvironnement global, qui dictera 1c choix 

5 des parametres en contrepartie d'une garantie d'une marge brute prd-dctcrmince (« (oiling 
agrccmc.nl »). Cette solution permet d*61imincr par mise en concurrence une dvcntucllc 
inefficacile d'un planiflcateur central, mais pose fondamentalemcnt 1c mcmc type de 
problfemes : les accords de ce type ne pcuvent jamais parfaitcmcnl prendre en comptc Ics 
contraintcs opdrationnelles locales. Elle conduit par ailleurs h communiquer & « Tassel 
10 manager » des informations confidentielles sur les contraintcs de fonctionnement du site. 

- Chaque producteur effectue son optimisation sur la base d'une moderation analytique de 
Penvironnement global et de son impact sur son 6quation de rendement, Ceci se fait au 
prix de fortes hypotheses simplificatrices qui, lorsqu'elles ne sont pas respectees, 
conduisent h des situations de crise (pertes tres elev^es pour le producteur, deficit de 

15 production globale, etc ...). 

- Chaquc producteur simulc tous les etats possibles de Tcnvironncmcnt (global ct 
specifique pour chacun des horizons de temps) corrcspondant aux sous-marchds 
d^electricitc. II calcule sa marge brute pour chaque £tat et chaque trajectoire dc 
parametres, et selectionne 1c vecteur de parametres optimisanl le couple « Esp6rance / 

20 VaR » sur sa marge brute. La simulation des atals possibles de Tcnvironncmcnt global 
impliquc la mise en oeuvre dc calculs complexes, mobilisant des moyens informatiques el 
des competences couteux, h F Evidence dupliques entrc les different* productcurs. 

La mise en oeuvre des enseignements de Tinvention permel & chaque producteur d'electricite 
dc choisir la troisteme approche, e'est-fc-dire la seule permcttant de faire une recherche 
25 precise de Toptimum local, tout en mutualisant la panic complcxc des calculs. Scion 
1'invention.: 

- II collccte rcpctitivcrncnt (avec une frequence corrcspondant h la frequence de eolation 
des prix sur le pool) laditc «smp s<5rie probabilisle » (W 4 f) muUi-va1u£c de description des 
pscudo-eiats probables futurs k sxt?ixp valeurs numeriques (W L f = A'/ global) qui lui est 

30 fournic par 1c centralise cmetteur (U). 

- Pour chaque « snip s£rie probabilisle », et pour chaque vecteur possible de parametres dc 
commandes, il tire une fonction aleatoire de defaut reflctant la probability de ddfaillance 
lechnique dc 1'installalion. U en deduil pour chaque 6tat el chaque vecteur dc paramctrc de 
commande une marge brute (A/B B (W//\X,)), calculi par une fonction arithmetique 

35 simple. 

- 11 en deduil par un tri simple classant pour chaque vecteur (K„) les tirages de sa marge 
brute (MB U ), unc fonction VaR„(K u ) ct unc fonction cspdrancc dc marge brute Emb u (k u ). 

Le probleme d'optimisation complexe en environnement aleatoire est ramen€ pour une 
centrale (S u iu s 1 ' ~ n) ) k un probteme d'optimisation simple: rechercher le maximum de 
40 EmbJKt) sous la contrainte VaR u (k tt ) < V max . 

Avantages de J'invention 

L 1 invention permet de minimiscr les moyens informatiques et humains n6cessaires h la mise 
en oeuvre de simulations probabilistes sur un ensemble de sites industriels distants soumis a 
un environncrneni aleatoire commun et vue d'aboutir nolammcnt ii : 
45 - unc optimisation de la production en respectant des contraintes reglemcntaircs, 

- une optimisation des couts de la production en controlant la survenue d'evenements 
defavorables, 

- un controle et une minimisation de la survenue d'evfenemcnts rarcs defavorables. 
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Un autre avantagc dc I'invention est de simplificr le calcul du risque global d'un groune de 
systemes industriels locaux distants. 

Un troisicme avantage de Pinvention est de permettre des caleuls dc simulation statistiquc 
locaux au niveau d'unc multitude de sites industriels a parlir dc tirages alcatoires mutualises 
5 sans quaucun des systemes locaux n'ait a transmetlrc d' informations au site centralis 
emetteur sur sa situation ou le r6glage de ses paramctrcs industriels. Ceci permct d'assurcr 
une meilleure confidentiality industriellc et commercial. 

En outre, Texcmplc duplication dans le domaine dc l'encrgie clectriquc donne ci-dessus 
permet de constater les avantages suivants pour un des producteurs industriels ■ 
10 - Le productcur peut se rapprocher de r optimum local, en mcttanl en ccuvre simplement ct 
avec peu de moyens locaux, les techniques d'optimisation les plus sophistiquces. Ceci se 
traduira par un temps dc retour plus rapide sur investi.ssemcnt ou par une baissc des prix 
au profit des consommateurs. 

II n'a pas besoin d'investir dans de couteux moyens de moderation et dc simulation de 
environnemcnt global. II peut se concenlrcr sur la simulation des paramelrcs specifiques 
locaux, independants de renvironnement global. 

II n'a pas besoin dc communiquer a un organisme central ses paramctrcs dc 
fonctionnemcnt. Dans un environnement concurrentiel, il peut done preserver la 
confidentialite dc sa strategic 



15 



20 



L'invention ayant maintenant etc decrite, et son interel just i fie sur des exemples detailles les 
demandeurs s'en reservent l'exclusivite pendant toute la duree du brevet, sans limitation autre 
que celle des lermes des icvendications ci-apres. 
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REVENDICATIONS 3 4 

1) Procede de simulation stochasdque repetitive liierarchisee el mulualisde par 
dchantillonnage de scenarios probables, pour la prediction probabilistc du comportement futur 
et/ou Poptimisation des paramfctres de rendement d'une multitude de syslimes industriels 
5 locaux distants (Sj u = M) ) : 

□ tres distanls les uns des autres, 

□ soumis * un environnement global commun (£) dvolutif, pouvant etre dccril 
numdriquement par p variables d'etat (XM (i=[ - p \ constiluant les composantcs d'un 
vecteur d'etat (X(r)) = (X,tl) (t = 1 ~ p \ 

10 □ ct dont refficacit6 de chaque systeme local (S u (u " l " n) ) peut en outre etre decrite par un 
indicateur numeriquc de rendement selon un modcle numdrique systemiquc 

local parametre (D u ) prdalablcmenl etabli (/„) = (D u (X u MO)\ dependant pour chacun 
des systemcs (S u ) et a chaque instant (r) : 

♦ de tout ou panic des variables d'etat (XKtf m * ^) du vectcur d^tat (X(/)), 
15 ♦ et d'un vectcur de q parametres (reglablcs) locaux d'aclion (X u ) = (A W |» /U); 
ce proc6d€ etant caractdrise en ce qu'en combinaison : 

□ au niveau d'un site centralist <*metteur (60, tres distant (de plusieurs kilometres, voire 
dc plusieurs centaincs de kilometres) d'une majoritd des n systcmes industriels locaux 
dislants (S/"" l " n) ) : 

20 4 on generc r£p£titivemem dc manicrc centralisee % des instants (r 0 ) successifs unc 
«smp serie de simulations probabilistes » (IV-/ <,3S| - •'">) mulli-vahide de 
description des variables (XAO iiu l " p) ) de renvironncment (£) constitute : 
a. soit de pseudo-6lats probables futurs des variables (X t ir) (its dc 
renvironncment (£"), 

25 b - soit dc pseudo-variations probables futures de ces variables } " p) ) dc 

Tenvironnement (£), 

communs a la multitude des /i (/; etant tres grand »1) systcmes industriels locaux 
distants (.?„"" ) decentralises, et simules ti un (on plusieurs ;//, avee m > I) 
horizons dc temps (//* = 1 "" r) > / c ) t 

30 ♦ ct on telctransmet rdpetitivement cettc « snip serie dc simulations 
probabilistes » (W,/) multi-valu£e de description des pseudo-dtats ou pseudo- 
variations probables futur(e)s des variables dc renvironncment, ce de manicrc 
rdpartie d^centralis^e : 
• du site cmetteur (£/), 
35 # en direction de V ensemble de la multitude des systcmes locaux distants (S„); 

u ct, de manicrc independante, au niveau de chacun des systemcs locaux \s u ) % 
majontairement tres distants du site centralise Cmetteur (U) t r<5p<5litivcmcn( : 

♦ on collectc ladite « snip s6rie probabiliste » (W,/) muhi-valuee de description des 
pscudo-eiats ou pscudo- variations probables futur(e)s, 

40 ♦ on cffcctuc un u-aitement num£rique statisiique local sp£cifiquc sur ladite s6rie 
probabiliste (W,/) re?ue pour predirc lc comportement futur probabiliste des 
systcmes industriels locaux (£„), 

♦ et (dventuellement) on regie le niveau du vecieur des paramitres locaux d % aclion 
<Xu) = iK\ K u ) de chaque systeme local (S v ) h une valcur optimalc calculee en 

45 function des criieres prectablis, afm noiammcnt d'aniicipcr le comportement 

probable futur et 1/cfficacite de chaque systeme (S u ) soumis h V environnement (E) 
au cours du ou des horizon(s) dc temps (t k ) retcnu(s). 

2) Procede scion la revendication 1 dc simulation stochastiquc repetitive hidrarchisee par 
cchantillonnagc de scenarios probables pour la prediction probabiliste du comportement futur 
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ci/ou I'optimisation des paramctrcs de rendement d'une muliiiude dc systemes induslriels 
locaux distants (S u ), ce precede ctant du type sptcifique consistant : 

□ au niveau d'un site expert (£/'), tres distant (dc plusicurs kilometres, voire de plusieurs 
centaines de kilometres) d'une majorilt des n systemes industriels locaux distants 

5 (S u {a ~ { " n) ), a faire prtalablement le choix : 

♦ des p variables d'environnement (X, (,s 1 ~ p) ) de 1'cnvironnemenl (E), 

♦ et d'un modele statisiique du comportcmcnl joint au temps (/) el par tranches de 
temps (Ar* ( * = l des variables d'etat (Xi(tf lr,] - p \ dtcrit par sa density de 
distribution jointe X p \ t\ Ar k )) : 

10 • justifies & partir d'observations historiqucs — il peut s'agir noiammeni de lois 

normales, log-normales, etc., ou plus generalement d'une distribution statistique 
jointe des p variables d'environnement (X,(t k ) 0 " L ~ p) ) ou de leurs variations 
tcmporcllcs (AX,(/ X ) (/= '"^ 

• el validt par des tests statistiqucs classiques adequats (moyennes, variances, 
15 correlations, sauts, comportemcnts extremes, ...), 

□ au niveau d'un site centralist tmetteur (LO, ties distant (dc plusicurs kilometres, voire 
dc plusieurs centaines de kilometres) d'une majorite des n systemes induslriels locaux 
distants (S u Ut=l " n) ): 

♦ & gencrer numtriquement repetitivemeni a des inslanis / c successifs, & partir du 
20 modele statistique retcnu (f^X\ X p \ t\ At k )) du comportcmcnl joint dans le temps 

(0 des variables d'etat (AT,), une «smp stric probabilisie » pj!sU ~"*) 

mutii-valuce de description des pseudo-etats ou pseudo-variations probables 
futur(e)s de renvironnement (E) constitute : 

• cas 1 : d'une « smp strie probabilisie » (Xij) dc pscudo-ttats formtc de sxm />- 
25 uplets d'cchantillons multi-values de pseudo-ctats futurs 

{Xiftd M "* % * X " M ""*)* comprcnant sxmxp valeurs name riques possibles des 
pscudo-ttats futurs probables des/? variables d'etat (A'0) f/ = 1 ' ct constitutes : 
• pour chacun des m (w > 1 ) horizons de temps Q k <k = * de la simulation, 

■ de .v (avec s ires grand » I) pseudo-etats futurs possibles (X fJ (td ijB "'Ydc 
30 chacunc des p variables d'etat (X ( (t k ) u = L , 

• cas_2 : et/ou tventuellement, d'une «.ww/> stric probabilistc varialionncllc » 
(AAT ( ,y ) formte de .vx/?7 /j-uplets d'cchantillons multi-valuts de pseudo-variations 
(AX tJ (t k ) ts:] - *■!«.»)) f u(urcs 5 partir dc la valeur courante des variables 
d'etat (AT/(/ 0 ) 0 = 1 ^ comprenant *xwx/> valeurs numtriaues de pseudo- 

35 variations futures probables des p variables d'tlat (Xj(f k ) lis { '~r) cl constituccs : 

■ pour chacun des m (r?i > I) horizons dc temps Q k iks de la simulation, 

■ dc 5 (avec s ires grand » 1) pseudo-variations numeriques futures possibles 
(A^<a) 0= L "* ) ) de chacune des /> variables d'etat (W" s| " /,) ), 

le procede ctant caracltrist en ce qu'en combinaison : 
40 □ au niveau du site centralist tmetteur (U) tres distant de la majorite des n systemes 
industriels locaux distants (S u (ua l ~ M) ), repetitivemeni a des instants (/ G ) : 

♦ on diffuse ladite «smp strie probabilisie » (W//) multi-valuee de description des 
pseudo-etats ou pseudo-variations probables futur(e)s, comprenant sxmxp pseudo- 
valeurs numtriques constitute par : 

45 a - < cas 0 ladite «smp strie probahilistc » (Xj) des pseudo-valeurs des variables 

d'etat (X i jCh)^ t ''+^ l "' s '* l '' M \ 
b. (cas 2) ei/ou ladite « smp stric probabilistc variationnclle » (AX tJ k ) des pseudo- 
variations (AX i j(!d ( ^ l ''*'^ u '' s '* mU ' JH \ 
ce de maniere dtcentraliste 
50 • du site emetteur (t/), 
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• cn direction de Tensemble de la multitude des systemes locaux dislams (S u ) ; 
□ et, de maniere independante, au niveau dc chacun des systemes locaux (S„), 
majoritairement tres distants du site centralist cmetteur ([/), rcpetitivemcnt : 

♦ on collecte laditc «smp seYie probabiliste » (Wj) multi-valuee de description des 
5 pseudo-6tats ou pseudo-variations probables futur(c)s h sxmxp valcurs numcriqucs 

♦ on choisit m'{m'< w) horizons de temps, rctenus parmi les dits horizons dc temps 

Ok ' " m )y 

♦ cvcntueliemenl, dans le cas (2) ou Ton reeoit uniquement ladite « smp seric 
10 probabiliste variationnelle » (AX,/) d'cchantillons variationncls (A*,//*)), on 

calcule Iocalement ladite «smp serie probabiliste » (#,/) des pscudo-dtats constitute 
de pseudo-cchantillons d'etais futurs ^ h i M "^ x '^ M "' m ^ correspondanie (au 
moins pour les m' horizons de temps (i k ) rctenus), 

♦ pour chacun des m' horizons de temps (/*) retcnus : 

15 # on apP'ique numeriquement Ie modclc systemique local (D u ) prealablernenl fixe 

afin de determiner dans chacun des «r scenarios j = 1,. ..,.y 1'indicateur de 
rendement : 

9 P° ur au moins une v aleur particuliere {\ u ) a = (A„ X« q ) du vecteur (K u ) des 

*v parametres locaux d'action, 

• et pour chacun des s p-uplets simulds de variables d'etat (X,./r*) (/= 1 " p \ 
correspondant audit horizon de temps (/*) rclcnu, de la seYie probabiliste de 
pseudo^tatC^./^) 02 '-^), 

et on determine les s pseudo-valeurs simulees dc Tindicalcur de rendement 
25 (UjOM) = (D^XjO^^ l -* correspondantcs, 

■ on rcgroupe et on ordonne les s valeurs correspondantcs dc rindicatcur dc 
rendement UujitkXuf* l "\ selon Ic niveau dc rendement (7^), 

" on caJ cule et on memorise Thistogramme (r w ) ainsi que hi distribution de 
30 probability (en nombre ou pourcentage) QY U (I U )) du niveau alleint dc Pindicaieur 

de rendement local (/«) dans la configuration (*„), selon 1c modele local (£>,,), 
• on en deduil par un calcul statistique la valeur d*au moins un parametre 
statisttque de rcglage choisi (V(X„)) (tel que : moyenne, ecarl type, Value-at- 
Risk, utilite, etc..) caracterislique : 
35 ' de la distribution de probability CV(K.f)) des niveaux de rendement (/„(/*X), 

■ cn fonction des variations du vecteur des parametres locaux diction 

(X U ) = (A W |,...,A W/ ), 

♦ (cvcntuellement) on effectuc Iocalement, au sein dc chaquc systemc local (S„), une 
optimisation numenque (selon des criteres predefinis) du (ou des) parametrc(s) 

40 stutistique(s) de rcglage choisi (Va*))> et on regie le niveau du vecteur des 

parametres locaux d'aclion a,) = (A„, X U(I ) de chaque systeme local (S u ) a unc 

valeur optimale calculee en fonction des criteres prcctablis, afin notamment 
d'anticiper le comportement probable futur et PcfficaciKS de chaque systemc (S u ) 
sourms h Tenvironnement (£) au cours du ou des horizon(s) de temps (/*) retcnu(s). 

45 3) Procede scion la revendication 2 de simulation slochastique repetitive hierarchisee par 
echantillonnagc de scenarios probables, caracterisc en ce qu'en combinaison : 
a au niveau du site expert (£/')» 

♦ on associc h chaque parametre local d'action (A„,/* = un code genenque (L,,) (ne 
dependant pas du systeme local (£,)), 
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♦ on dtablit pour chaque code (L A ), et on actualise regulicrement, lin «6valuatcur de 
rendemcnt » (n,,), c'est-a-dire un programme (materialise sous forme d'un fichier 
numerise) permettant de calculer, en fonction du p-uplel des variables d'etat 
(X) = (X,'- 1 et du parametre local d'action (X u ,,) y la composantc correspondante 

5 du rendement (7 M/f ) = (Tl^A^X)), 

♦ dc telle sorte que, pour chaque systfcme local (5„), ledil modfele systcmique local 
(D u ) consiste simplement a additionner les composantes (I uh ) de l'indicateur de 
rendement correspondant a chaque param&tre local d'action (X uh ) : 

0.tt..x)»£n,tf»,x) (*) 

10 □ au niveau d'un site centralist tmetteur (U): 

♦ on g£n6rc numdriquement et on emet, repctitivement & des instants (r 0 ) successife, h 

partir du modele statistique retenu {fJX J( p \ /; At)) du compoilcment joint dans le 

temps (0 des variables d'etat (AT/), la «smp scric probabilistc » (W,/) multi-valu£e 
de description des pseudo-dtats ou pseudo- variations probables futur(c)s > 

15 ♦ et on dmct, a Papparilion de chaque nouvcau code (L h ) ou dc chaque modification 
de r<$valuateur dc rendemcnt (U,,\ la paire (L/M fl ) sous forme d'un fichier 
num6ris6, 

□ au niveau dc chaque systfcme industriel local distant (S lt ) : 

♦ on calculc : 

20 # P° ur c ^acun des m' {m'< m) horizons de temps (/*) dc la simulation relenus h 

r echelon local 

• pour chacun des / = ! s (avee s trcs grand » i) pseudo-etats numdriques 

possibles (Xj -* ) ) recus, ou reconsiitues h partir dc la s£ric des pseudo- 
variations (AXij). 

25 # P our chacun des h = 1 q param&tres locaux d'action (A*,,), 

• la composantc de rendement (/„/,/) = (n^A^X/)), 

• puis, par sornmation, selon la formulc (*), Pindicateur dc rendemcnt (/„/), 

• le parametre statistique retenu (V(k„)), 

♦ on applique le vecteur dc paramctres locaux diction (X.J qui optimise la valcur du 
30 parametre statistique (V). 

4) Precede selon la revendication 2 dc simulation slochastique repetitive hieYarchisce par 
echantillonnage dc scenarios probables, caracterise en cc qu'en combinaison : 
a au niveau du site expert (U') % 

♦ on associc h chaque parametre local d'aciion Uj h = 1 "«>) un code gdntrique (L f( ) (nc 
35 dependant pas du systeme local (S u )) t 

♦ on ctablit pour chaque code (£,„), et on actualise rdgulierement, un «<5valuateur de 
rendement marginal » (ir h ), c'est-a-dire un programme (materialise sous forme d'un 
fichier numcrisd) permettant de calculer, en fonction du /duplet des variables d'etat 
(X) = (A'/ / 1 la composante correspondante du rendement marginal (/'i,) = 

40 (atfX)), 

♦ de telle sorte que, pour chaque systftmc local ledit modele systcmique local 
(D u ) consiste simplement a effectuer un produil seal aire cntre Ic vecteur de 
parametre locaux d'action (X u ) = ( Xj 1 s » "*>) et le vecteur des rendements 
marginaux (I*) : 

45 DM = (k u \n (**) 

□ au niveau d'un site centralise dmettcur (U) : 

4 on gSnerc numdriquement et on emet, repctitivement a des instants (t 0 ) successifs, a 

partir du modele statistique retenu {f^X } X p \ /; A/))du comportemcnt joint dans 1c 
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temps (f) des variables d'etat (X,), la «swp sdric probabiliste » (W,f) multi-valu£e 
de description des pscudo-6tats ou pseudo-variations probables futur(e)s, 

♦ et on emet, a P apparition de chaque nouveau code (L h ) ou de chaque modification 
de I'evaluateur de rendement (*,,), la paire (L h> n h ) sous forme d'un fichier numerisS, 

5 ♦ et on calcule : 

• pour chacun des in horizons de temps (/*), 

• pour chacun desy= l,..„s (avec s iris grand » 1) pseudo-^tats num<5riques 
possibles (X/ 0=, "*V 

• pour chacun des = 1 param&tres locaux action (A,,,,), 
10 • la composante de rendement marginal (/ '„/) = (*/><X/)), 

♦ et on diffuse la « smq serie probabiliste » (W' h j) multi-valuee de description des 
pscudo-dtats (/'/,/) ou pseudo-variations probables futur(e)s '(fiJ'rf) , 

□ au niveau de chaque systfcme Industrie] local distant (S u ) : 

♦ on collecte la « smq sdrie probabiliste » (W',,j) multi-valu6c de description des 
15 pscudo-^tats (7*/).ou pseudo-variations probables futur(e)s (A/*/) 

♦ puis, par sommation, scion la formulc (**), on cn deduit la sdric des indicatcurs de 
rendement (/„/) et le paramfetre statistique rctcnu (V(X it )) y 

♦ on calcule, puis on applique le vecteur de paramfctres locaux d'action (X u ) qui 
optimise la valeur du paramfctre statistique (V). 

20 5) ?voc6d6 scion la revendication 2 de simulation stochastiquc rdpdlitivc hidrarchisdc par 
echantillonnagc dc scenarios probables, caractdrisd en ce qu'en combinaison : 

□ au niveau d'un site centralist £metteur (U): 

♦ on genfcre numdriqucmcnt rdp6titivement a des instants (/ c ) successifs, h partir du 
module statistique rctcnu (f<AXi % .„,X p \ t\ Af)) du comportement joint dans 1c temps 

25 (/) des variables d'etat (X, ( ' s '* une serie probabiliste variationncllc (AX,/) 

formcc dc sxm /duplets d'cchantillons multi-values dc pseudo-variations 
(AX"///*) 1 " 8 ■■■f , '/-i---*.* B «---*rt) d'etats futurs h partir dc T<Stat courant des variables d'dtat 
(X / </ 0 ) <f ~ , " ,>) ), comprenant sxmxp valeurs numcriqucs dc pseudo-variations futures 
probables des p variables d'etat (Xfaf* l " p) ) et constitutes: 

30 • pour chacun des /// (m > 1) horizons de temps - m> ) de la simulation, 

• dc s (avec s tres grand » I) pseudo-variations numdriqucs possibles 
(AX,-//*)*' = des p variables d'dtat l " p \ 

♦ on diffuse de maniere repetitive ladite «smp strie probabiliste variationnetle » 
(AY,/) des pseudo-variations d'6tals (AXjj(uf° l - p ' Jsl - 5 ' M ~ M \ ce de maniere 

35 decentralisee : 

• du site emetteur (U), 

• en direction de Tensemble de la multitude des systemes locaux dislants ($„), 

□ et, dc manifcrc independante, au niveau dc chacun des systemes locaux (S u ) : 

♦ on collecte et on memorise de maniere repetitive h chaque instant (r 0 ), scion unc 
40 frequence determinee, voire & chaque modification de Tune des variables (X,), la 

valeur courante des variables d'etat (X l {i 0 f ml - py ), 

♦ on collecte ladite «smp s<Srie probabiliste variationncllc » (AX,/) des variations 
d'dtals (&XiA**) (M "+ M '' J * M '''*>) en fonction de Pincrdmcnt de temps Ar k , emise 
par le site central (U), 

45 ♦ on calcule localemenl ladite «smp siric probabiliste » (X,/) de pseudo-etats 
(Xiflk) M "+ m ^ 1 ''*^ 1 "^) aux instants futurs (/*) = (/ 0 + A/jt), corrcspondant a la serie 
variationncllc par calcul numdrique incremental du type : 

ou toute autre formule (addition, exponenticllc, etc. . .) faisant intervenir : 
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• ladite « smp serie probabiliste variationncllc » (AX,/) emise par lc site centralist 
dmetteur (LO, 

• et la valeur courante des variables d'dtal (X,</) (,= l " p) ). 

6) Procddd selon la revendication 2 de simulation stochastique repetitive hidrarchisee par 
5 echantillonnage de scenarios probables, caractdrisd en ce qu'en combinaison : 
□ au niveau du site centralist dmetteur (£/): 

♦ on gdndre numdriquement rdpdtitivement a des instants (Q successifs, a partir du 
modclc statistiquc rctcnu (faPC u ...J( p \ /; AO) du comportcment joint dans le temps 
(/) des variables d'etat (X,), une « smp sdric probabiliste variationnelle » (AX//) 
formde de sxni 77-uplets d'dchantillons multi-values dc pseudo-variations 
(AX^/*/ ,B, "^ J " l "* ,rBl -"" ,) ) d'etats futurs a partir dc Tctat courant des variables d'dtat 
(X^/J 0 "* 1 "^), comprenant sxmxp valeurs numcriques de pseudo-variations futures 
probables des p variables d'etat (X t it k f= lp) ) ct constitutes : 

• pour chacun des m (m > 1 ) horizons de temps (// te! - w) ) de la simulation, 

15 • de s (avee s tvbs grand » I) pseudo- variations numdriqucs possibles 

(AX u 0 k f= '■••*) des p variables d'etat (X.W* l ' p \ 

♦ on collccte de manidre repetitive ct on memorise la valeur de Tetat courant des 
variables d'etat (X ( -(/ 0 )'' = ~ p \ 

♦ on diffuse dc manidre rdpdtitive et dcccntralisee du site dmetteur (LO cn direction dc 
20 Pcnscmble de la multitude des systemes locaux distants (S„), a la fois : 

• la «smp sdric probabiliste variationncllc » (AX,/) des pseudo-variations 
(AX^/,)^ 1 -^^ 1 -^ 81 --)), 

• et la valeur courante des variables d'etat (X,(f 0 )" = ' ' :P \ 

7) Pvoc6d6 selon la revendication 2 de simulation stochastique rdpdtitive hidrarchisee par 
25 echantillonnage de scenarios probables caractdrisd en ce qu'eri combinaison : 

a prcalablement, et en difftrents points de renvironnemem (E), on installe une sdrie de r 
capteurs - r) ) observatcurs de l'dtal d'environncment (/T), cmcttant a chaquc 

instant (/), r mesures dc Tetat de Tenvironncmcnt (p u O) 0t " 1 * r) ), 

au niveau du site conccptcur «/')> on definit prcalablement un modclc numerique 
desenptif d'environncment (<D) f ddfinissant Tcnsemble dc p variables d'etat 
(X/(0" =1 " p) ) = (0(^i(/),...,iU f (/))) en fonction des mesures (JJt h (J) ih ^ x "' r) ) (avee 
gdndratement p <> /•), en sorte que l'ensemblc des variables d'environnement 
(X,(/) - - ) decrivent physiquemcnl avee une precision acceptable, 5 chaque instant 
d'dchantillonnage (/*), la globalitd de l'dtat dc renvironnement (£)', 
35 □ au niveau du site centralisd dmetteur (60: 

♦ on collectc a distance ct on mdmorise les r mesures d'environnement (/i/,(/) f * * 1 * r) ), 
ces /• mesures constituant a chaque instant d'dchantillonnage (/*) les composantes 
d'un vecteur temporel de mesure d'environncment (M(t k )) = (jj,,(t k f' = , " r) ) 1 

♦ el on calciile, a chaque periode d'dchantillonnage (f k ) successive, a Taidc du modele 
40 numdnquc descriptif d'environnement (fl>) et du vecteur de mesures (M(;)) les p 

variables d'etat (XM Uml "+ } ) = (O(p,(/),...,M0)). 
8) Procddd selon la revendication 2 de simulation stochastique rdpdtitive hidrarchisde par 
echantillonnage de scenarios probables, caractdrisd en ce qu'en combinaison : 
a au niveau du site centralise dmetteur (U): 
45 ♦ on maintient et on met $ jour un historique des valeurs des p variables d'dtal 
(X,(/) ' ' " p ) dc Tcnvironncment global commun (£) au cours du temps. 

♦ on gdndre numdriquement rdpdtitivement a des instants (r 0 ) successifs, a partir : 
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• du modele statisiique retcnu (f^X,,...JC p ; t; At)) du comportenient joint dans le 
temps (/) des variables d'etat (X," = 

• el de eel historique mdmorise des valeurs des p variables d'ttat (X,{t)° = 1 ■''>) 

la «sm P serie probabiliste » (W,f) multi-valuce de description des pseudo-etats ou 
5 pseudo-variations probables futur(e)s de 1'cnvironncment (£), 

♦ on diffuse ]adkc« snip s6rie probabiliste » (W//) multi-valuee de description des 
pseudo-etats ou pseudo-variations probables futur(c)s de manierc deccntralisee du 
site cmetteur (60 en direction de I 'ensemble de la multitude des systfcmes locaux 
distants (S u ). 

10 9) Procedc scion la rcvendication 8 de simulation slochastique repetitive hierarchisec par 
echantillonnage de scenarios probables, caractcrise en ce qu'en combinaison ■ 
□ au niveau du site centralist dmetteur (U) : 

* fv/T-V 1 -*)" 1 Ct on met * J° ur un historique des valeurs des p variables d'etat 
(X,(0 ) de 1'environnement global commun (£) au cours du temps, 

15 ♦ on generc mimcriqucment repetitivemenl a des instants (/„) successifs, a parlir : 

• du modele statisiique retcnu (fc(* If ... Jf,; t; At)) du comportemcni joint dans le 
temps (/) des variables d'etat (Xf"° 

• ct de cct historique memorise des valeurs des p variables d'etat (X,it)" = '-'") 

la «smp serie probabiliste » (W,j) multi-valuce de description des pseudo-etats ou 
-o pseudo- variations probables fulur(e)s de 1'environnement (£), scion la methode ditc de 

« Monte-Carlo », par tirages pseudo-aleatoires des pseudo-valeurs (W,h grScc h un 
generateur de nombres situe sur le site centralise emctteur (U), voire sur le site 
conceptcur (£/'), reposant sur unc mdthode : 

• soit « pscudo-aleatoire » de nature ergodiquc, 

25 * soU « " discrepance faible » (suites de Sobol, Hamcrsley, etc.), 

• soil sur touic autre methode produisant des nombres asympiotiquement repartis 
selon une distribution specifiec (par cxemplc uni forme sur un inlervallc), 

♦ on diffuse ladite « stnp serie probabiliste »> (W-J) multi-valuce de description des 

pseudo-etats ou pseudo-variations probables futur(c)s de manierc deccntralisee du 

M sue emettcur (U) en direction de Penscmblc de la multitude des systemes locaux 

distants (£„). " " 

10)Procede scion la revendication 8 de simulation slochastique repetitive hierarchisec par 
echantillonnage de scenarios probables, caractcrise en ce qu'en combinaison : 
□ au niveau du site centralist emetteur (U): 

♦ ??,"<?l" li Tj\. on n,el h -» our un historique des valeurs des p variables d'etat 
(A,W ) de 1 cnvironnement global commun (E) au cours du temps, 

♦ on genere numeriquement repetitivemenl a des inslants (/„) succcssifs, a' parlir • 

• du modele statisiique rctenu (a*, .,*„; ,; At)) du comportcment joint dans le 
temps (0 des variables d'ttat 1 '"), 

• et de cet historique memorise des valeurs des p variables d'etat (X,itf = 1 ^) 

la «smp serie probabiliste » (W,f) multi-valuce de description des pseudo-etats ou 
pseudo-variations probables futur(e)s de renvironncment (FT), scion la methode 
r w r< "! S T qUe * d ' cs,imatjon d 'eiats futurs probables des p variables d'etat 
[5 X,) ' } recons,ruit s * Partir de Petal courant (X,(t„f- 1p) ) et de variations 
45 h.stonqucs passees (AX^^^) pour un ensemble de dates passees 

( ' •' ' < '«>)• ce en affeclanl h chaque etat reconstruit une probability de realisation 

pouvant dependre de rage de la variation historique correspondanle, 

♦ et on diffuse ladite «smp serie probabiliste » (W,/) multi-valuce de description des 
pseudo-etats ou pseudo-variations probables fuiur(e)s dc manierc decentralisee du 
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site dmetteur (60 en direction de Tenscmble de la multitude des systfemcs locaux 
distants (S tt ). 

11) Procddc selon la revindication 8 de simulation stochastiquc repetitive hidrarchisdc par 
dchanlillonnage de scenarios probables, caracldrisd en ce qu'en combinaison : 

5 □ au niveau du site centralist dmetteur (U) : 

♦ on maintient et on met k jour un historique des valeurs des p variables d'etat 
(XiU) 0=l ' p) ) de l'environncment global comtnun (£) au cours du temps, 

♦ on gdn&rc numdriquement repdtitivement h des instants (/ 0 ) successifs la «smp serie 
probabiliste » (IV//) multi-valuee de description des pseudo-dlats ou pseudo- 

10 variations probables futur(e)s de l'environncment (E), h partir : 

• du module statistique retenu <fa(X u ...J( p \ t; A/)) du comportement joint dans le 
temps (/) des variables d'etat (Xj), 

• et de cet historique mdmorisd des valeurs des p variables d'etat (Xi(t) (iss 

• mais en introduisant volontairement dans la « smp sdrie probabiliste » (W-/) des 
15 hiais de comportement statistique joint des variables d'etat (XK0 iiu 1 " p) ) ' vis-ft- 

vis dc leur comportement statistique historique, 
les series biaisdes reproduisant ainsi des « scenarios catastrophes », afin de tenir 
compte de V impact desdits scenarios catastrophes sur les sysicmcs distants (SJ, 

♦ on diffuse ladite « smp sdrie probabiliste » (Wj) multi-valude biaisdc de description 
20 des pseudo-dtats ou pseudo-variations probables futur(e)s dc manifere decentralisdc 

du site dmetteur (U) en direction de Tensemble dc la multitude des sysicmcs locaux 
distants (S u ). 

12) Procedd scion la revindication 2 de simulation stochastiquc rdpdtitive hidrarchisee par 
cchantillonnage de scenarios probables, caracterisd en ce qu'en combinaison : 

25 u certains des systfemes distants (S u ) envoient au systcme central (U) des requetes sous 
forme de parametres de scenarios catastrophes scion des biais specifiques definis par 
cux, 

□ au niveau du syslfcme central : 

♦ on gendrc en fonction de ces requetes, evenluellemcni en fusionnant ccllcs qui sont 
30 idemiques, des «smp series probabilistes » (W,/) multi-valudcs de description des 

pseudo-etats ou pseudo-variations probables futut(e)s, biaisdes en fonction des 
parametres specifies par les systfcmes distants, 

♦ et on envoic du site dmetteur (U) vers chaque sysl&me local distant (S u ) ladite « smp 
sdrie probabiliste » (W,f) multi-valuee biaisee de description des pseudo-etats ou 

35 pseudo-variations probables futur(e)s qui correspond a sa requcte. 

13;Procddc selon la revendication 2 de simulation stochastiquc repetitive hidrarchisee par 
cchantillonnage dc scenarios probables, caracldrisd en ce qu'en combinaison : 

□ chaquc systdme distant (S u ) alterc lui-meme tout ou panie des «smp sdrics 
probabilistes » (W t f) multi-valuees de description des pseudo-dlats ou pseudo- 

40 variations probables futur(e)s, de fa?on & cn faire des series biaisdes reproduisant des 
scenarios catastrophes dont les parametres sont : 

♦ soit spdeifids au niveau du systcme central (If) et diffuses en direction des systdmcs 
locaux (S u )y 

♦ soil definis au niveau des systemes locaux (S u ). 

45 14)Procdde scion la revendication 2 de simulation stochaslique rdpdtitive hidrarchisdc par 
cchantillonnage de scenarios probables, caracterisd en ce qu'en combinaison : 

□ au niveau du site^ centralisd dmetteur ({/) on gdndre numdriquement la « smp sdric 
probabiliste » (Wif) multi-valuee dc description des pseudo-etats ou pseudo-variations 
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probables futur(c)s de Tenvironnement (£) dc manitrc pdriodique, selon une frequence 
d'echantillonnage (<p). 

15) Proctde selon la revendicalion 2 de simulation slochastique repetitive higrarchiste par 
echantillonnage de scenarios probables* caracttrist cn ce qu'en combinaison : 

5 □ au niveau du site centralist tmetteur (U): 

o on diffuse ladite «smp strie probabiliste >> (W,/) multi-valuee de description des 
pseudo-ttats ou pseudo-variations probables futur(c)s de manifcre ddcenlraliste du 
site tmetteur (U) cn direction de I'ensemble de la multitude des systtmes locaux 
distants (S u ) par Emission sous une forme « broadcast », nolamment par Emission 
10 radio, 

□ et, de manitre independante, au niveau dc chacun des systcmcs locaux (5 tt ), 
majoritairement trfes distants du site centralist emetteur (£/), rtpttilivement : 

o on collecte ladite «smp strie probabiliste » (Wif) multi-valute de description des 
pseudo-ttats ou pseudo-variations probables futur(e)s & sxnixp valeurs num^riqucs 
15 (Wif = Xif ou AXif). par reception quasimcnt cn continu dc cctle Emission 

« broadcast » en provenance du site centralise emetteur (U). 

16) Procede selon la revendication 2 de simulation stochastiquc rtptlitive hidrarchiscc par 
echantillonnage de scenarios probables, caracltrise en ce qu'en combinaison : 

□ au niveau du site centralist tmetteur (U): 

20 o on diffuse ladite «smp strie probabiliste » (Wif) multi-valuee dc description des 
pseudo-ttats ou pseudo-variations probables futur(e)s, de mantere dtcentnilisee du 
site tmetteur (U) en direction de P ensemble de la multitude des syslfcmcs locaux 
distants (S„) par Emission sous une forme dite « PP aval », c*cst-<Vdirc par Emission 
sur un rescau (R) (notamment de type Internet) du site centralist tmetteur (U) 5 
25 Tadresse sp6cifique de chacun des systcmcs industricLs locaux distants (S u {u * '"" n) ), 

a et, de manifere indtpendante, au niveau de chacun des systcmcs locaux (5 U ). 
majoritairement trfes distants du site centralise emetteur (U). rtpctilivcment : 
<> on collecte, ladite «smp strie probabiliste » (Wif) multi-valute de description des 
pscudo-etats ou pseudo-variations probables futur(e)s h sxmxp valeurs numdriques 
30 (Wif = X,f on AX,f) t par enregistremcnl succcssifs dans la mtmoire (//„) de 

Tordinatcur local (C u ). 

17) Procdde selon la revendication 2 de simulation stochastique repetitive hicrarchisec par 
echantillonnage dc scenarios probables, caractdris6 cn cc qu'en combinaison : 

□ au niveau du site centralist tmetteur (U): 

35 0 on diffuse ladite « snip strie probabiliste » (Wif) multi-valuee dc description des 
pseudo-ttats ou pseudo-variations probables futur(e)s dc manidrc dtcentraliste du 
site tmetteur (U) en direction de P ensemble de la multitude des systemes locaux 
distants (S„) par misc & disposition de ladite «smp strie probabiliste » (Wif) multi- 
valute sur un serveur (2) disponible 
40 o sur un rSseau (R) (notamment de type Internet), 

o par accts « PP » par chacun des systtmes industriels locaux distants (Sj" = 1 , 
u et, dc manure independante, au niveau dc chacun des systfcmcs locaux (S u ) 
rtpttitivemenl : 

o on collecte ladite « smp strie probabiliste » (Wif) multi-valuee dc description des 
45 pseudo-ttats ou pseudo-variations probables futur(e)s h sxtnxp valeurs numeriques 

(W, f j = Xq ou AX,/), par accts sous unc forme dite «PP amont », e'est-a-dire par 
acefcs via 1c rtscau (R) (notamment dc type Internet) au site centralist tmetteur (I/) h 
son adresse specifique, par acquisition par transferl (notamment de type « FTP ») de 



39/48 



BNSOOC1D: <FR 2B08909A 1.. L > 



H 0 2808909 
ladite «smp sErie probabiliste » (VV,/) multi-valuEe ct par cnrcgislrement dans la 
memoire (H u ) de Tordinateur local (CJ. 

18) ProcEde selon la revendication 2 de simulation stochastique r^p^titive hiErarchisec par 
echantillonnage dc scenarios probables, caractErise en ce qu'cn combinaison : 

5 □ an niveau du site centralist Emetteur (U) y on gEnEre numEriquement rcpeiitivement h 
des instants (r 0 ) successifs, ii partir du module stalistique retenu (fcPCu...JC P \ t; A/)) du 
comportement joint dans le temps (/) des variables d'etat (X*' unc «smp sErie 

probabilistc » (Wi/) mulli-valuEe de description des pseudo-Etats ou pseudo-variations 
probables futur(e)s de Tenvironnement (E) constitute de $ (avec s trcs grand » 1) 
10 pseudo-valcurs numEriques des p composantes (S LJ (! k f =] ' ' ,,Jsl ' s) ) pour chacun des m 
(m> 1) horizons de temps (//*" , *" m> ) l sous la forme d'un fichier dc sxmxp donnEes 
alpha-numEriques successives subdivisE : 

♦ en m « sous-fichiers temporels » successifs de niveau 1 regroupant les pseudo- 
valeurs d'un memc horizon de temps 0 k ik = l " M '\ 
15 ♦ chaque « sous-fichier tempore 1 » incluant une succession dc s « sous-fichiers de 
tiragc » de niveau 2 constitutes chacun d'un /duplet de pseudo-valeurs. 

19) Proccdc selon la revendication 2 de simulation stochastique repetitive hierarchisee par 
echantillonnage de scenarios probables, caracterise cn ce qu'en combinaison : 

□ au niveau du site centralise emetteur (£/), on gEnfcre numEriquement rdpeti live men t & 
20 des instants (/ 0 ) successifs, une «smp sErie probabiliste » (W,/) mulli-valucc dc 

description des pseudo-Etats ou pseudo-variations probables futur(c)s dc 
Tcnvironncment (E) en separant les m « sous-fichiers temporels » successifs dc niveau 
1 regroupant les pseudo-valeurs d'un memc horizon de temps (t k ikw 1 •' H> ) par un 
indicateurd*horizon de temps (rj fc w ; = 1 ""' ) ) corrcspondant, 
25 □ et, au niveau de chacun des systemes locaux (£„), on distingue les uns des autres les 
differents dits sous-fichiers temporels d'une meme « snip sErie probabilistc » (W,./) 
multi-valuEe collectee en idemifiant leurs indicateurs specifiqucs d v horizon de temps 

20) Procedc selon la revendication 2 de simulation stochastique rEpclilivc hierarchisee par 
30 Echantillonnage dc scenarios probables, caractErisE cn cc qu'en combinaison : 

□ au niveau du site centralist Emetteur (U), on genfcre numEriquement repEtitivement h 
des instants (/ 0 ) successifs, unc « snip serie probabiliste » (W,f) mulli-valuee de 
description des pseudo-etats ou pseudo-variations probables futur(c)s de 
Tenvironnemeiu (£). cn sEparant les s « sous-fichiers de tirage » successifs dc niveau 2 

35 (d'un memc fichier temporel) regroupant chacun un />-uplet de pseudo-valcurs par un 
indicaieur de tirage (%f Jm '— J) ) correspondant, 

□ et, au niveau de chacun des systemes locaux (S u ) % majoritairement trEs distants du site 
centralise Emetteur (U) % on distingue les uns des autres les differents « sous-fichiers de 
tirage » successifs dc niveau 2 (d'un mSmc fichier temporel) regroupant chacun un p- 

40 uplei de pseudo-valeurs d'une memc « snip sErie probabilistc » (IV,/) multi-valuce 
collcctcc cn identifiant leur indicateur spEcifique de tirage (rf?~ 1 , ' ) ). 

21) Proccdc selon la revendication 2 dc simulation stochastique repetitive hidrarchisee par 
echantillonnage dc scenarios probables, caractErisE en ce qu'en combinaison : 

u au niveau du site centralise emetteur (£/): 
45 ♦ on diffuse dans sa globalite ladite « snip serie probabiliste » (W { f) multi-valuee de 
description des pseudo-Etats ou pseudo-variations probables futur(e)s de maniEre 
dEcentralisEc du site Emetteur (60 en direction de Tcnscmble de la multitude des 
systemes locaux distants (S„) par Emission sous une forme « broadcast », 
notammcnt par Emission radio, accessible a tous les systEmcs locaux (S tl ), 
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♦ et, dc manifcrc independante, au niveau dc chacun des systemes locaux (S u ) t on 
installe et on utilise un fihre stlectif prtrtgle dc mise & disposition d*unc fraction 
(WijY preddfinic (notamment acquise par un abonncmcnt specifique) de ladite 
«smp serie probabiliste » (W,/) multi-valute de description des pseudo-etats ou 

5 pseudo-variations probables futur(e)s, colleclee par reception en provenance du site 

centralist emetteur (60 • 

22) Proctdt scion la revendication 2 de simulation stochastique repetitive hitrarchiste par 
tchantillonnage de scenarios probables, caracttrist en cc qu'en combinaison : 

□ au niveau du site centralise tmctteur (60> on diffuse dc mantere selective (selon une 
10 r£gle d'attribution prtdtfinic) : 

♦ divcrscs fractions differences ((Hfy*)/*' l " Jrt ) prtdefinies (notamment acquises par 
un abonnement sptcifique) dc ladite «smp strie probabiliste » (Wjf) multi-valucc 
de description des pseudo-etats ou pseudo-variations probables fulur(e)s, 

♦ de manidrc dccentraliste du site emetteur (60 en direction dc plusieurs groupes 
15 (^ (A = 1 "'*') attributifs dc la multitude des systemes locaux distants (S u \ chaque 

groupe (G g ) recevant une fraction particuliere ((W,/)^ 

♦ par Emission sous unc forme dite « PP aval >\ e'est-ft-dire par Emission sur un rtseau 
(R) (notamment dc type Internet) du site centralist tmetteur (60 a 1'adresse 
specifique de chacun des systemes induslricls locaux distants (S,/" = 

20 □ et, de manitre indtpendante, au niveau de chacun des groupes (G u { * = 1 " M ) attributifs de 
systemes locaux (S tt ) % majoritairement trts distants du site centralise emetteur (60, 
rtpetitivement : 

♦ on collccte ladite fraction prtdtfinie ((W,/)*) dc la «smp serie probabiliste » (W,/) 
multi-valuee dc description des pseudo-ttals ou pseudo-variations probables 

25 futur(c)s par enregistrements successifs dans la mtmoire (//„) dc 1'ordinateur local 

(C). 

23) Proccde scion la revendication 2 de simulation stochastique repetitive hitrarchiste par 
cchantillonnagc de scenarios probables, caracttrist cn cc qu'en combinaison : 

□ au niveau du site centralist tmetteur (60 : 

30 ♦ on diffuse ladite «smp strie probabiliste » (VV,/ ) multi-valuee dc description des 
pscudo-etals ou pseudo-variations probables fu(ur(c)s dc manifere dccentraliste du 
site emetteur (60 en direction de V ensemble dc la multitude des systemes locaux 
distants (S„) par mise a disposition de ladite «smp strie probabiliste » (W,./) multi- 
valuee sur un scrvcur (Z) disponible sur un rtseau (R) (notamment dc type Internet) 

35 P ar accts PP par chacun des systemes industriels locaux distants (S t/ (u = 

♦ cn mcltant dc manitre selective & disposition de dif ft rents groupes (G v ( * s 1 " M ) dc 
systemes locaux (S u ) sur le servcur (Z) diverses fractions Iimitatives ((VV,/) V ) 
prtdefinies (notamment acquises par un abonnement sptcifique) de ladite «smp 
serie probabiliste » (IV,-/) multi-valuee de description des pseudo-ttats probables 

40 futurs, ce notamment par Tintermtdiairc d'un systtme d'idcnlificateurs et dc mots 

de passe transmis par chaque systeme local (S u (u = 1 * rt) ), 

□ et, de manitre independante, au niveau de chacun des systemes locaux (S u ) y 
rtpttitivement : 

♦ on collccte ladite fraction .((W,/),) prtdtfinic attribute dc ladite « snip strie 
45 probabiliste » (W,/) muki-valutc de description des pscudo-ttats ou pseudo- 
variations probables futur(e)s par acces sous unc forme dite « PP amont », c'est-<V 
dire par accts via le rtseau (R) (notamment de type Internet) du site centralist 
emetteur (60 & son adresse sptcifique, acquisition par transfer! dc ladite fraction 
((Wq\) prtdtfinie attribute dc ladite « smp serie probabiliste » (W,/) multi- 

50 valute et par enregistrement dans la mdmoire (H u ) de 1'ordinateur local (C M ). 
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24) Proceed selon la revendicaiion 2 de simulation stochastique repetitive hi6rarchis<5c par 
echamillonnage de scenarios probables, caractcrisc cn ce qu'en combinaison : 

□ on choisit prealablement n'(n'< n) systemes industries locaux, retenus parmi les n dits 
systemes indusiricls locaux (S u (u = 1 " n \ 

5 a on choisit prealablement m' (m'<m) horizons de temps, retenus parmi les m dits 
horizons de temps 0k ks l " l \ identiques pour tous les systemes locaux, 

□ au niveau de chacun des n' systemes locaux (S„) retenus, pour chacun des w' horizons 
de temps retenus : 

♦ on applique numeriquement le modele syslcmiquc local (DJ prealablement fixe afin 
10 dc determiner dans chacun des s scenarios 0 = 1, • Tindicateur de rendement 

(/„/ = D U (X U% (Xjj(t k )) : 

• pour au moins une valeur particulicrc (X u ) du vecieur des paramelres locaux 
d'action, 

• et pour chacun des s /duplets simules de variables d'etat (X t j(i k ) (i = l '~ p \ 
15 correspondant audit horizon de lemps (/<) retenu, de la serie probabiliste de 

pseudo-etats QCjihf* 1 ~ s \ 

• et on determine les s pscudo-valcurs simuiees dc Tindicateur de rendement 
(Iuj(hX)) = (D u (\ u ,XjO k ))V= 1 - 5) )correspondantcs, 

♦ on emct la «serie de simulations locales » constiluec des s pseudo-valeurs simuiees 
20 dc Tindicateur de rendement (I UJ <t k X u )) = (D it (X u Xj0 k )f = '"^ corrcspondantes du 

systeme local de maniere centralisante en direction d'un site centralisateur (£), 

□ au niveau du site centralisateur (£), 

♦ on collecte les sxn' pseudo-valeurs simuiees d'indicateurs dc rendement 
ihjihXu)^' , "* > ) en provenance des chacun des n' systemes locaux (S H ) cl pour 

25 chacune des s simulations, 

♦ on calcule, pour chaque simulation (/= l,...,.v) un indicatcur de rendement global 
(J } ) de Tensemblc des systemes locaux a partir d'un modele devaluation global 
(/,) = (W^X)' 1 ^-"^). 

♦ on regroupe et on ordonne les s valcurs corrcspondantes de Tindicateur dc 
30 rendement (/,(/* ^u...,*.,,-) 0 ' 51 [ '\ selon son niveau, 

♦ on calcule et on memorise Thistogramme (D ainsi que la distribution dc probabilite 
(cn nombrc ou pourcentage) 0P(7)) du niveau atleint de Tindicateur de rendement 
global (J) dans la configuration choisie des vecteurs de paramelres locaux d'action 
(A) = (X tt {u=:l M \ selon les modeles locaux (/)„) et le modele devaluation global (D) 

35 retenus, 

♦ on cn deduit par un calcul statistiquc la valeur d'au moins un paramctrc statistique 
dc rdglagc global choisi (V(A)) (tel que : moyenne, ccart type, Value at Risk, utilitd, 
etc..) caracteristique : 

• dc la distribution dc probabilite (xCku-.-X, J)) des niveaux de rendement 
40 gIobaux(7(^;X, \ n )), 

• en fonction des valcurs du multi-vecieur (A) = (k\ % ...Jk n ) des vecteurs dc 
parametres locaux d'action (k tt ) y 

♦ (cvcntuellement) on effectue globalcmcni unc optimisation numeriquc (selon des 
criteres predefinis) du parametre statistiquc de reglage global choisi (V(A)) par 

45 rapport au multi-vecteur (A) = (k \, t ) cl on transmel a chaque systeme local (S„) 

le niveau correspondant dc son vecteur optimise des parametres locaux d'action (k u ) % 
afin notamment d'anticiper globalement lc comportcmcnt probable futur et 
Tefficacite de Tenscmble des systeme locaux (S u ) soumis a Tcnyironnement (£) au 
cours du (ou des) horizon(s) de temps (/*) retenu(s). 
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25)Proc6d6 selon la revendication 2 de simulation slochastiquc rdp^titivc hidrarchisdc par 
echaniillonnage de pseudo-6tats futurs probables de variables de marche, pour la prediction 
probabiliste du comportement fulur et/ou r optimisation des paramfctrcs de rendement d'une 
multitude de n (n etant tits grand »1) d'opiratcurs financiers locaux (S„ (w = U n) ) : 
5 □ decentralises et tres distants les uns des autrcs, 

□ poss6dant chacun un portefeuille (f u (W5C l "" > ) : 

♦ constitu6 d'actifs financiers (0/, (A s u 

♦ donl la valeur (I J" " ] ~ n) ) ddpend des mouvements du marchc des capitaux (E), 

♦ et dont la valeur (I u {u c ] ~ n) ) peut etre calculec par combinaison : 

10 • d'un vccleur dc variables dc marchd (X(/)) = (X t {tf * 1 " p) ) communes variant au 

cours du temps (/), & partir duquel il possible, pour chaque opdratcur financier 
local distant (S t{ ) x dc reconstituer avee une precision acceptable la valeur de tous 
les actifs financiers pouvant €tre presents dans son portefeuille (P„), scion un 
module local devaluation (D„), (ces variables dc march£ pouvant repr6senter 

15 des prix d'actifs ou de commodit&s, des taux d*int£ret ou des coefficients de 

courbe de taux, des volatililes ou des coefficients dc surfaces dc volatility des 
indicateurs de defaut d'un 6metteur ou d'unc contrcparlie, ou tout autre nombre 
apportant unc information pertinente sur Petal des marches de capitaux), 

• ct, pour chaque operateur financier local distant (S tt )> d'un veclcur de quantitds 
20 (k tt ) = (X** 1 *" '"•«*) de chacun des actifs financiers (ft ( * s , ~* ) ) detenus dans son 

portefeuille (PJ, par exemplc scion la combinaison lindaire : 

U0 = £A u/ ,/r^X(0A) (***) 

ou (n u (uX{t) y e it )) rcprdscntc lc prix dc Vacuf(e„) a Hristant (0 selon le modelc 
local (D„) lorsque lc march6 est dans Tdtat X(/), ou une autre formule pouvant 
25 fairc apparaTtre une non-lin6aritd par rapport aux quantity (A,,/,), pour tenir 

comptc notammcnt du risque de liquidite; 
ce proedde etant du type consistant : 

□ au niveau d*un site expert (t/')» trfcs distant (de plusieurs kilometres, voire de plusicurs 
centaines de kilometres) d'une majoritc des n opdratcurs financiers locaux distants 

30 (S M { " U 1 >d faire prcalablcment le choix : 

♦ des p variables dc marche ((Xif = 1 ~ p) ) decrivant Tenscmblc du marchc (E), 

♦ et d'un modele statistique du comportement joint au temps Q) et par tranches de 
temps (A//* = 1 " p,) ) des variables de marchd m (XK0 iim l " p \ d<5crit par sa densite dc 
distribution jointe (J^X { X p ; t; A/J) : 

35 • justifie h partir d'observations historiqucs — il peut s'agir notammcnt dc lois 

normalcs, log-normales, etc., ou plus gdncralement d'une distribution statistique 
jointe des p variables d'environncmcnt (A',</<) (f = , * v>) ) ou dc leurs variations 
temporelles (&XAi k ? m ] - p \ 

• et valide par des tests statistiques classiques adequats (moyennes, variances, 
40 correlations, sauts, comportements extremes, ...). 

□ au niveau d'un site centralise dmetteur (£/)» tres distant (dc plusicurs kilometres, voire 
de plusieurs centaines de kilometres) d'une majority des n systemes induslrieis locaux 
distants (S u ( * =, "" ) ) : 

♦ h gdndrcr num£riquement r6p6titivemcnt a des instants (/ 0 ) successifs, h partir du 
modelc statistique retenu (f<JtXu...X P \ t\ A/*)) du comportement joint dans le temps 
(/) des variables dc marche " unc « snip s6rie probabiliste » (W { f) multi- 
valuee de description des pseudo-etats ou pseudo-variations probables futur(c)s du 
marche (£) constitute : 
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• cas 1 : d'une « smp strie probabilisle » (X,f) de pseudo-etats formee de sxm p- 
uplets d'echantillons multi-values de pseudo-etats futurs 
{X i Jt£ < ^ l ''*'fc } " J '** l '- M \ comprenant sxmxf) valeurs numtriques possibles des p 
variables de marche (XjOf* i " p) ) ct constitutes : 

5 • pour chacun des m (m > 1) horizons de temps (f t ( *° dc )a simulation, 

• de s (avec s trfcs grand » 1; pscudo-ttats futurs possibles des /; variables dc 
marche '"*V 

• cas 2 : et/ou eventuellement, d'une «smp strie probabiliste variationnelle » 
(AXjj) formte de sxm p-uplets d'echantillons multi-values de pseudo-variations 

10 (AXi//*) 1 '" 1 "^^ 1 "**" 1 -^) futures a partir de 1'etat courant des variables de 

marche (X&o)*' * 1 "' p \ comprenant sxmxp valeurs numtriques de pseudo- 
variations futures probables des p variables d'etat (JWtt) 0 " - 1 et constitutes : 

• pour chacun des m (m > 1) horizons de temps (i k {k = 1 " m) ) de la simulation, 

■ de s (avec s trts grand » 1) pseudo-variations numtriques possibles 
15 (AXi/ik)^ 1 '"*) des p variables de marche QCfaf* L 

Ic proctdc ttanl caracttrist en ce qu'en combinaison : 

□ au niveau du site centralise emetteur (U) tr6s distant dc la majoritc des n optratcurs 
financiers locaux distants (S u iu * L repetitivement a des instants (/ 0 ) : 

♦ on diffuse ladite « snip strie probabiliste » (VV,/) multi-valuce dc description des 
20 pseudo-ttats ou pseudo variations probables futur(c)s, comprenant sxmxp pscudo- 

valcurs numtriques, constitute par : 

a. (cas 1) ladite «s?np strie probabiliste » (Xtj) des pseudo-valeurs des variables 
dc marche (Xi/h) 1 **' ■^ , ■> B '■■■ , •* ,a, ■■ m \ 

b. (cas 2) el/ou ladite « snip strie probabiliste variationnelle » (AXjf) des pseudo- 
25 variaUohs(AXi i/ (f A ) tfBl "^> l "**" L - J " ) ) 1 

cc dc manifcrc dtccntralistc 

• du site emetteur (£/)» 

• en direction de Tensemble de la multitude des operateurs locaux distants (S h ) ; 

□ et. de mantere indtpendante, au niveau de chacun des operateurs locaux (S u ). 
30 majoritaircmenl trfcs distants du site centralist tmeltcur (I/), rtpelitivement : 

♦ on collcctc ladite « smp scric probabiliste » (W,j) multi-valucc dc description des 
pscudo-ttats ou pseudo-variations probables futur(e)s a sxmxp valcurs numtriques 

♦ on choisit m'Qn'< m) horizons dc temps, rctcnus parmi les /// dits horizons de temps 
35 (h sk : l ' m \ 

♦ eventuellement, dans 1c cas (2) ou Ton recoil uniqucment ladite « snip strie 
probabiliste variationnelle » (AX,/) d'echantillons variationnels (A*,-//*)), on 
calculc localcmcnt ladite « smp scric probabiliste » (Xif) des pseudo-etats constitute 
de pseudo-tchaniillons d'ttats futurs (XiJt k ) tM '-+^ l ''* M '" my ) corrcspondanlc, (au 

40 moins pour les w' horizons de temps (/*) rctcnus), 

♦ pour chacun des w' horizons dc temps (r k ) retenus : 

• on calcule numeriquemenl scion le modclc local devaluation (£>„) 
prtalablement fixt, dans chacun des s sctnarios (J = l»...,.v), la valeur globalc 
(/*/') = (D u (\ { „Xj(t k )) du portefeuille (P u ) : 

45 ' P°w au moins une valeur pariiculitre (X w )= (A„ t ,...,A„,) du vecteur (K u ) des 

quantites d'aclifs detenus, 
■ et pour chacun des s p-uplets simules de variables d'ttat (JO//*)), 
corrcspondant audit horizon de temps 0 k ) retcnu, dc la strie probabiliste dc 
pseudo-etat^/x) 0 " 1 ~ s \ 
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et on determine les s pseudo-valeurs simuldes du ponefeuille 
= DuiKiXjitkfi* 1 - ^correspondantcs, 

■ on regroupe et on ordonne les s valcurs corrcspondantcs (/„,//* A w )° = 1 " s) ) du 
portefeuille (P u ) y 

5 «on calcule et on memorise la distribution de probability (en nombre ou 

pourccntagc) V ? w (/ M ) dc la valcur (/„) du portefeuille local (P u ) dans la 
configuration (?v„), scion le modele local (D„), 
• on en deduil par un calcul statistiquc la valeur d'au moins un param&trc 
statistique de reglage choisi (V(\ u )) (tel que: moyennc, dcart type, Value-al- 
10 Risk, utilite, etc ...) caracteristique : 

• de la distribution dc probability W((k u )>0 des valcurs O u Q k% \ u ), 

■ ct du vecteur (X b ) = (A„ K g ) ddcrivant lc portefeuille (P tt ), 

♦ (eventuellcmcnt) on effectue localement, au niveau de chaquc op6rateur local (S u ) t 
une optimisation numdrique (selon des entires prdddfinis) du (ou des) paramfetrc(s) 

15 statistiquc(s) de reglage choisi (V(\ u )) t ct on ajustc la composition du poficfcuille 

(K) = (/Ui, ...Xt U ) de chaque operateur local (S u ) a une valeur optimale calculdc en 
fonction des critics prcetablis, afin notammcnt d*anticipcr Ic comportcmcnt 
probable futur ct le rendement de chaque portefeuille (P u ) soumis au marchd (E) au 
cours du ou des horizon(s) de temps (/*) retenu(s). 

20 26)Procddc scion la revindication 2 de simulation stochastique repetitive hidrarchisce par 
cchantillonnagc de pseudo-etats futurs probables dc variables de marche, pour la prediction 
probabiliste du comportcmcnt futur et/ou 1' optimisation des paramdtres de rendement d % unc 
multitude dc n (n dtant trfes grand » 1) d'opdratcurs financiers Iocaux (5 M (WBt , - #rt ) : 
□ ddccntralisds et tres distants les uns des autres, 

25 □ possddant chacun un portefeuille (Pj um ] "' n) ) : 

♦ constitud d'actifs financiers (0 h <hs 

♦ dont la valeur unilaire (/*** = depend des mouvements du marchd des capitaux 
(E), et est une fonction d'un vecteur de variables de marche (X(/)) = (X/(0° e l "* > ) 
communes variant au cours du temps 0), & partir duqucl il possible, avec une 

30 prdcision acceptable, dc rcconsliluer la valeur de tous les actifs financiers pouvanl 

eire presents dans les portefeuilles (P u ) des actcurs locaux, selon un modcie commun 
devaluation (D), (ces variables de marchd pouvant rcprcscntcr des prix d'actifs ou 
de commodates, des taux d'interel ou des coefficients de courbe de taux, des 
volatilites ou des coefficients de surfaces de volatility des indicaieurs dc defaut d'un 

35 emetteur ou d'une contrepartie, ou tout autre nombre apportant une information 

pertinente sur Pdtat des marches de capitaux), 

♦ et dont la valeur (I u iu s 1 •" ) ) peut dtre calculee par produil scalaire entre : 

• le vecteur des valeurs unitaires (/') '= (/V* * 

• et, pour chaque operateur financier local distant (S w ), d'un vecteur de quantitds 
40 (K) = (Kh ih ~ '"■**) de chacun des actifs financiers (fi h iHm 1 "« } ) ddtenus dans son 

portefeuille (P u ) : 

D^X) = (k it \ /') (***) 
oO (/'/,) = (7r/,(/,X(0,ft,)) reprdscnte le prix dc ractif(0/,) a Tinstant (0 selon le 
modele commun (D) lorsque le marchd est dans 1'dtat (X(r)), 
45 ce proeddd dtant du type consistant : 

a au niveau d'un site expert (£/')> trcs distant (de plusicurs kilometres, voire de plusieurs 
centaines de kilometres) d'une majoritc des n operateurs financiers locaux distants 
(S* " M ), a faire prdalablement le choix : 

♦ des p variables de marche ] - p) ) ddcrivant Pcnscmblc du marche (ZT), 
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♦ et d'un modele statistique du comportemcnt joint au temps (/) et par tranches de 
temps (At k a = 1 •'">) des variables de marche (X,<if m ] " p \ dccrit par sa densitd de 
distribution jointe (f ( £X\ X p \ t\ A/*)) : 

• justify h partir d'observations historiqucs — il pcut s'agir nolammcnt dc lois 
5 normales, log-normales, etc., ou plus gendralemcm d'unc distribution statistique 

jointe des p variables d'environnement (Xj(t k ) (i = ou dc leurs variations 
temporalis (AftC/*)""? ~+\ 

• et validd par des tests statistiqucs classiqucs adequats (moyennes, variances, 
correlations, sauts, comportemenls extremes, etc.), 

10 □ au niveau d'un site centralist dmettcur (10, tres distant (de plusicurs kilometres, voire 
dc plusieurs centaines de kilometres) d'une majorite des n systdmcs industriels locaux 
distants (Sj" m ' ~ n) ) : 

♦ h gdndrcr numdriquement rdpdtitivcment h des instants </ 0 ) successifs, h partir du 
modele statistique rctcnu (fdXu : *.<X P \ K A/a,)) du comportcment joint dans le temps 

15 (/) des variables de marche (Xf* * i " />) ), une « snip serie probabilistc » (W,/) mulli- 

valuee dc description des pseudo-dtats ou pseudo-variations probables futur(c)s du 
marche (£), constitute : 

• cas n°l : d'une «s?7jp sdrie probata lisle » (X}/) dc pscudo-etats formcc de sxm 
/7-uplcts d'echantillons multi-values dc pseudo-etats futurs 

20 com p renant s xnixp vaicurs numcriques rcprdscntant : 

■ pour chacun des m (m > 1) horizons de temps (i* tt " I -"' ) ) dc la simulation, 

■ s (avec 5 ties grand » 1) pseudo-diats futurs possibles des p variables de 
marche .(X,W" 

• cas n°2 : d'une «smp sdrie probabiliste variationncllc » (AAT,/) dc pseudo-etats 
25 formdc dc sxni /duplets d'cchantillons multi-values dc pseudo-variations futures 

(AXij(! k ) {i=i - PJxl ' s,k ^-'' m) ) i comprenant sxmxp valeurs numdriques rcprdscntant : 

■ pour chacun des w (m > 1) horizons de temps (/< a = 1 de la simulation, 

■ ,v (avec s trds grand » J) pseudo-variations futures possibles des p variables 
de marche (XKh) u w 1 * 

30 ♦ & calculer pour chaque vecteur (X,/) de pscudo-etal des variables de marche unc 

« smq sdrie probabiliste » (W/,/) multi-valucc dc description : 

• cas n°l : des pseudo-prix probables futurs d'actif 

• cas n°2 : et/ou des pseudo-variations (A/*//*)'*" 1 ""'* /rL *** " w> ) qui en 
35 dccoulent 

le proctde etant caracterist en ce qu'en combinaison : 

□ au niveau du site centralise tmctteur (if) irts distant de la majorite des n optrateurs 
financiers locaux distants (S u {u = , " Jrt ). rdpetitivement il des instants (t 0 ) : 

♦ on diffuse, ladite «smq sdrie probabiliste » (W'hj) multi-valude de description des 
40 pseudo-etats ou pseudo variations probables futur(e)s, comprenant sxmxq pseudo- 

valeurs numdriques constitute par : 

• (cas 1) ladite « smq serie probabiliste » (// v ) des pseudo-valeurs des prix des 
actifs (/' A / (fcssK "*>■■ 

• (cas 2) et/ou ladite « smq sdrie probabiliste variationncllc » des pscudo- 
45 variations (A/*/**" 1, > n >\ 

ce dc manifcre decentralisee 

• du site emetteur (60, 

• en direction de Tensemble de la multitude des operateurs locaux distants (S u ) ; 

□ el, de manitre independante, au niveau de chacun des operateurs locaux (S u )> 
50 majoritairement tres distants du site centralise dmettcur (£/), rdpetitivement : 
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0 on collccie ladite « smq s£rie probabiliste » (TV'/,/) multi-valuce de description des 
pseudo-6tals ou pseudo-variations probables fulur(c)s h sxmxq valcurs numeriques 

(vn/ = /;./ouA/;/), 

o on choisit m'(m'< m) horizons de temps, rctenus parmi les m dits horizons dc temps 
5 Ok ' 

o dventuellement, dans le cas (2) ou Ton recoil uniquernent ladite « smq s<Srie 
probabiliste variationnelle » (A/;/) d^chantillons varialionncls (A/*//*)), on 
calcule Iocalemcnt ladite «smq s£rie probabiliste » (/',,/) des pseudo-Stats 
constitute de pseudo-6chantillons ductals futurs U%J0k) ihsX ' •* s -* =, -" rt ) 
10 correspondantc, (au moins pour les w' horizons de temps (/*) rctenus), 

o pour chacun des m ' horizons de temps (/<) rctenus : 

o on calcule numdriquement scion le modele local devaluation (D„) 
prealablement fixe dans chacun des s scenarios (/ = U-..,s) la valcur globale 
(/«/) = (K\ /;*) du portefeuille (P u ) : 
15 D P our ™ moins unc valcur particulicrc du vecteur (X u ) = (A^i,.:.^) des 

quantitds d'aclifs detenus, 
° et pour chacun des s* duplets simulcs dc prix d'actifs (// v {/*)), correspondant 
audit horizon de temps (t k ) retenu, dc la s6ric probabiliste dc pseudo-dtats 

20 ct on determine les s pseudo-valeurs simul<5es du portefeuille (J uJ 

0* ;A w )) = (A, J r/)^' ^) correspondantcs, 

o pour la valcur (X u ) = (A wl A u „) du vecteur des paramfctrcs locaux d'action 

rctenus, 

° on rcgroupe et on ordonne les s valeurs corrcspondanles (I u ,{t k Xif s l "' s) ) du 
25 portefeuille (P„), 

D on calcule et on memorise Fhistogrammc (r„) ainsi que la distribution dc 
probability (cn nombre ou pourcentage) QV u (F tt )) dc la valeur (/„) du portefeuille 
local (P u ) dans la configuration (Xu), selon lc modfclc local (£>„), 
o on cn ddduii par un calcul statistiquc la valeur d'au moins un paramctre 
30 statistiquc dc reglagc choisi (V(X tl )) (tel que : moyenne, dcaii type, Value-at- 

Risk, utilite, etc. . .) caractcristiquc : 
° de la distribution de probability 0F„(X, 4 ,/ M )) des valeurs (J u (t k XX 
° en fonction du vecteur (X u ) = (A wl , ...,A W</ ) decrivant le portefeuille (P w ). 
❖ (evcntuellcmcnt) on effectue Iocalemcnt, au niveau de chaque optrateur local (5 W ), 
35 une optimisation numtrique (selon des criteres prfdefinis) du (ou des) parametre(s) 

statistique(s) dc rdglage choisi (V(X U )\ ct on ajustc la composition du portefeuille 
(K) = (A*|, de chaque opdrateur local (S u ) & une valeur optimale calculee en 

fonction des criteres preetablis, afin notamment d'anticiper le comportement 
probable futur et lc rendemcnt de chaque portefeuille (P„) soumis au marche (£) au 
40 cours du ou des horizon(s) de temps (/*) retenu(s). 

27)Proc<§de selon la revendicalion 25 de simulation stochastique repetitive hierarchisee par 
echantillonnagc de pseudo-6tats futurs probables de variables dc march€, pour la prediction 
probabiliste du comportement futur et/ou ['optimisation des parametrcs de rendement d'une 
multitude de n (n etant tres grand » 1) d*op<5rateurs financiers locaux , ce proc<*d(5 etant 
45 caracteris£ en ce que : 

□ pour cvaluer la valeur (/„/) du portefeuille (/>„) d'un opdrateur financier (S u ) h la date 
future (/*) selon ie module (D u ) et dans T&at du marche simulc (X//*)) correspondant 
au scenario (/), on prend en compte dans le modele (£>„) les &ats de marche simules 
dans le meme scenario (/), aux dates ulldrieures (i/'"** 1 - -0 ), selon le principe suivant : 

47/48 
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♦ on dgfinit pour chaque actif financier (d h (h " 1 une quantity maximale (\> h ) qu'ii est 
possible de (raiter (5 Taehat ou & la venle) dans une journde, 

♦ si aucune des quantites (\ uH ) d'actifs (6 h ) detenus dans les portefeuilles (P u ) n'cxcfcdc 
en valeur absolue la valeur limite (v,,), alors on applique la formule (***) pour ^valuer, 

5 & la date (/*) et dans la simulation (j) la valeur (/„/) du portefeuillc (P u ) : 

♦ si une des quantity (X uh ) d'actifs (G h ) ddlenus dans un des portefeuilles (P u ) cxcfcde en 
valeur absolue la valeur limite (v h ) t alors on remplace la quanlite (A uh ) par (w,), 
affectde de son signe, a la date (t k ) et on reporte le reliquat (X uh ± v h ) sur la date 

10 suivante (/*+,) en utilisant le vecteur de variables de marche (Xj(h+\)) ; 

♦ si cette quanlitd 6tait encore en valeur absolue supdrieure h la limite <Y„), on la fixe b la 
valeur maximale (v, f ), avee le mtae signe, et on reporte & nouveau le reliquat (X uh ± 
2\> h ) sur la dale suivante (t M ) en utilisant le vecteur de variables de marchd (X/0 k+2 )), 
et ainsi de suite... ; ■ 

15 ceci, notamment pour tenir compte du risque de liquiditd imposant des contraintcs dc borne 
supdricurc & la quantite dc chaque actif (0,,) que Ton est en droit de valoriser selon le 
modele local (Z>„) h une date donnee (/*), et entramant un report dventuel de la position 
vinucllc rcsuliante dans chaque simulation sur des dates ulldricures. 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 9 
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Les revendications presentes contiennent 
des options , des expressions vagues et des 
termes conditioned. El les manquent 
depressions techniques, Elles sont 
dirig€es d'une maniere confuse et 
trompeuse. Un manque de clarte s'ensuit, 
De plus il est impossible de determiner 
I'objet pour laquelle une protection est 
demandee. Le manque de clarte, et 
1 'impossibility de determination sont 
telles qu'une recherche significative de 
I'objet des revendications devient 
impossible. Par consequent, aucun rapport 
de recherche ne peut etre etabli pour la 
prosente demande. 

Seulement pour exemple non limitatif, les 
revendications 1 et 2 contiennent , entre 
autres choses, les options, expressions 
vagues ou termes conditioned suivants: 



- tres distants 

- un environement global commun evolutif, 
pouvant etre decrit num§riquement 

- peut en outre etre decrite par un 
indicateur 

- on g£nere . de maniere central is^e 

- une «smp serie 

- pseudo-etats probables futurs 

- pseudo-variations probables futures 

- horizons de temps 

- un traitment num§rique statistique local 
spScifique 

- un traitement ... pour pr§dire le 
comportement futur probabiliste 

- et (eventuel lement) on regie le niveau 
... a une valeur optimale calculee en 

-/-- 
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